Ground-state phase diagram of the one-dimensional half-filled extended Hubbard model 
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We revisit the ground-state phase diagram of the one-dimensional half-filled extended Hubbard model with 
on-site (U) and nearest-neighbor (V) repulsive interactions. In the first half of the paper, using the weak- 
coupling renormalization-group approach (g-ology) including second-order corrections to the coupling con- 
stants, we show that bond-charge-density-wave (BCDW) phase exists for U ~ 2V in between charge-density- 
wave (CDW) and spin-density-wave (SDW) phases. We find that the umklapp scattering of parallel-spin 
electrons disfavors the BCDW state and leads to a bicritical point where the CDW-BCDW and SDW-BCDW 
continuous-transition lines merge into the CDW-SDW first-order transition line. In the second half of the paper, 
we investigate the phase diagram of the extended Hubbard model with either additional staggered site potential 
A or bond alternation 5. Although the alternating site potential A strongly favors the CDW state (that is, a band 
insulator), the BCDW state is not destroyed completely and occupies a finite region in the phase diagram. Our 
result is a natural generalization of the work by Fabrizio, Gogolin, and Nersesyan [Phys. Rev. Lett. 83, 2014 
(1999)], who predicted the existence of a spontaneously dimerized insulating state between a band insulator and 
a Mott insulator in the phase diagram of the ionic Hubbard model. The bond alternation 5 destroys the SDW 
state and changes it into the BCDW state (or Peierls insulating state). As a result the phase diagram of the 
model with S contains only a single critical line separating the Peierls insulator phase and the CDW phase. The 
addition of A or 5 changes the universality class of the CDW-BCDW transition from the Gaussian transition 
into the Ising transition. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.10.Pm, 71.30.+h 



I. INTRODUCTION 



It is well known that a one-dimensional (ID) spin sys- 
tem has instability to dimerization that changes the system 
into a nonmagnetic insulating state, the so-called spin-Peierls 
state.' Indeed the spin-Peierls state is realized in many sys- 
tems including quasi-one-dimensional organic compounds^-^ 
and the inorganic material'* CuGeOs, and its properties have 
been studied extensively both experimentally and theoreti- 
cally. Of particular interest is a situation in which a dimerized 
state appears spontaneously due to strong correlations and 
frustration.^ A well-known example is the frustrated spin-i 
Heisenberg chain with nearest-neighbor, Ji, and next-nearest- 
neighbor, J2, antiferromagnetic exchange interactions, where 
a spontaneously dimerized phase is realized for J2 > J2c — 

0. 24J1.'' Other systems of cuiTent interest are quasi-one- 
dimensional electron systems in organic materials, where the 
spin-Peierls state appears due to strong electron correlation at 
half fiiiingMi2iiaiiiii^iiiii and at quarter filling 

Recently it was pointed out by Nakamura' and co-workers 
that a spontaneously dimerized state occupies a finite param- 
eter space in the ground-state phase diagram of the ID half- 
filled Hubbard model with the nearest-neighbor repulsion V, 

1. e., the extended Hubbard model (EHM). This spin-Peierls 
state is often called bond-charge-density-wave (BCDW) state 
or bond-ordered- wave state. The appearance of the BCDW 
state in the purely electronic model is nontrivial and has at- 
tracted much attention from theoretical point of view. To ap- 
preciate this surprising result, let us consider some limiting 
cases. In the limit of weak nearest-neighbor repulsion V, or in 
the half-filled Hubbard model with only the on-site Coulomb 



repulsion U, the ground state is in the Mott insulating state 
where the spin sector exhibits quasi-long-range order of spin- 
density wave (SDW); we call it the SDW state. In the opposite 
limit of strong V, the ground state of the half-filled EHM has a 
long-range order of the charge-density wave (CDW); we call 
this state the CDW state. Furthermore, in the atomic limit 
where the electron hopping t is ignored, the CDW state ap- 
pears for U < 2V whereas the uniform state coiTesponding to 
the SDW state is stable for U > 2V in one dimension. Strong- 
coupling perturbation theory in t has established that a first- 
order phase transition between the SDW state and the CDW 
state occurs at [/ ~ 2y.'^-'^-^"-2' As for the weak-coupling 
regime, perturbative renormalization-group (RG) approach or 
g-ology led to a similar conclusion that the ground state at half 
filling is either in the SDW state or in the CDW state with a 
continuous phase-transition line at U — 2^.'* Thus, it had 
been considered for a long time that the ground-state phase 
diagram of the EHM at half filling has only two phases, the 
SDW and CDW states, and that the order of the phase tran- 
sition at J7 ~ 2V changes from continuous to first order at a 
tri critical point which was speculated to exist in the intermedi- 
ate coupling regime This common view was revised 
by the Nakamura's discovery that the BCDW state exists at 
U ~ 2V in between the SDW and CDW phases in the weak- 
coupling regionfii which is supported by recent large-scale 
Monte Carlo calculations. ^^-^^ Related studies of the dimer- 
ized state in the E HM with additi onal correlation effects can 
be found in Refs. I27l28l29l3(l3ir . 

A related and still controversial issue of cuiTent interest 
is whether or not a spontaneously dimerized phase exists in 
the ID Hubbard model with alternating site potential, the 

so-called ionic Hubbard modeI.3^.33,34,35,36,37,38,39,40,41,42,43,4V5 
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This system was introduced as a simple minimal model for the 
neutral-ionic transitions observed in quasi-one-dimensional 
organic materials^^i^iffi and for ferroelectric perovskites42iSi 
Obviously the model has two insulating phases. The ground 
state is (i) a band insulator with the CDW order when the stag- 
gered site potential is much larger than the on-site repulsion or 
(ii) a Mott insulator with quasi-long-range SDW order when 
the staggered site potential is negligible. Early exact diagonal- 
ization studie£'*2i^2ii of small systems have found a transition 
between the two phases and also reported dramatic enhance- 
ment of the electron-lattice interaction by strong electron cor- 
relation near a boundary between the band insulating phase 
(the BI state) and the Mott insulating phase (the SDW state). 
Mostly through bosonization analysis of the ionic Hubbard 
model, Fabrizio, Gogolin, and Nersesyan recently argued^ 
that a phase of a spontaneously dimerized insulator (SDI) in- 
tervenes between the ionic insulating phase (band insulator) 
and the Mott insulating phase. The SDI state is closely re- 
lated to the BCDW state mentioned above. Earlier numerical 
studies^i^^i^^i^2i^2ii have drawn contradictory conclusions as 
to whether the SDI phase exists or not, but more recent nu- 
merical studies find two phase transitions and the SDI phase 
in betweeniSiiSii^ Nevertheless there still remain unresolved 
issues on the critical properties near the quantum phase tran- 
sitions. 

In this paper we give supporting theoretical arguments for 
the existence of the spontaneously dimerized insulating states 
in the ID half-filled extended Hubbard model with and with- 
out staggered potentials. We adopt the standard bosonization 
approach and perform both perturbative RG analysis valid in 
the weak-coupling regime and semiclassical analysis which 
is expected to give a qualitatively correct picture even in the 
strong-coupling regime. This paper is organized as follows. 
Sections|n]and|ni]are devoted to the analysis of the standard 
EHM, i.e., the system without the staggered potential. Some 
of the results of this part are already presented in Ref. 52. In 
Sec. HH we introduce the model and reformulate the weak- 
coupling theory, the y-ology, to include higher-order correc- 
tions to coupling constants. We bosonize low-energy effec- 
tive Hamiltonian and derive the renormalization-group equa- 
tions. In Sec. Hn] we determine the ground-state phase dia- 
gram. First, from the perturbative RG analysis we show that 
the BCDW phase occupies a finite region near the U = 2V 
line in the weak-coupling limit. Next, from the semiclassi- 
cal analysis we argue that the umklapp scattering of parallel- 
spin electrons destabilizes the BCDW phase and gives rise to 
a bicritical point where the CDW-BCDW and SDW-BCDW 
continuous-transition lines merge into the CDW-SDW first- 
order transition line. Finally, combining the perturbative RG 
equations with the semiclassical analysis, we obtain the global 
phase diagram of the ID EHM. In Sec. llVl we study the ID 
EHM with the staggered site potential. We take the same strat- 
egy as in the previous sections and perform a semiclassical 
analysis of the bosonized Hamiltonian. With the help of the 
perturbative RG analysis we obtain the global phase diagram 
that indeed has the SDI phase. We find that the BCDW phase 
of the EHM is continuously deformed to the SDI phase upon 
introducing the alternating site potential. In Sec.|V] we study 



the ID EHM with additional bond dimerization, but without 
the staggered potential. This model exhibits a quantum phase 
transition between a dimerized Peierls insulator and a CDW 
state. Section lVTl is devoted to conclusions, and details of the 
technical calculations are given in Appendixes. 



II. EXTENDED HUBBARD MODEL 

In the first half of this paper (Sees. HI] and IIIH . we con- 
sider the standard ID EHM which has on-site, U, and nearest- 
neighbor, V, interactions. The Hamiltonian is given by 

+ U^n^,]nj,i+V^njnj+i, (2.1) 



where nj,a- = Cj.o-S'.'t ~ h = + "jM' ^l.a 
notes the creation operator of an electron with spin cr (= |, 
I) on the jth site. We assume repulsive interactions, i.e., the 
coupling constants U and V are positive. Note that the Hamil- 
tonian has global SU(2) spin symmetry. Following the previ- 
ous studies on models with correlated-hopping interactions,^^- 
we consider the CDW, SDW, BCDW, and bond-spin-density- 
wave (BSDW) phases as potential ordered ground states at 
half filling. They are characterized by the order parameters 

CcDW = (-!)■' {rij.i + Uj^i), (2.2a) 

OsDW = (-1)^ (ri,,T - (2.2b) 

Obcdw = (-l)^(4^c^.+i^^ + c] ^c^.+i^^ + H.c), (2.2c) 

Obsdw = (-l)^(4|C^.+i^i- - cl^c^+i,i + H.c). (2.2d) 

The order parameter of the BCDW state corresponds to the 
Peierls dimerization operator We note that the BCDW 
state can be also regarded as the p-density-wave state, as 
the order parameter of the BCDW state can be written as 

EjCbcdw oc Efc,^sin(A;a)4,<,c^+(^/^)^^, where Ck,a = 
^-1/2 ^-ikRj^,^^ ^jjjj j^, _ j-^. jjjg jattice spacing, 
N: the number of sites). The BSDW state describes a site-off- 
diagonal SDW stateiS 



A. p-ology approach 

The hopping t generates the energy band with dispersion 

£fe = — 2t cosfca, where the Fermi points are at fc = ztkp — 
±7r/2a at half filling. In order to analyze the low-energy 
physics near the Fermi points, we introduce a momentum 
cutoff A(0 < A < kp) and divide the momentum space 
into the three sectors (Fig. ^ (i) fc e R, (ii) k £ L, and 
(iii) k {RU L), where R = [kp - A, kp + A] and 
L = [—kp — A,—kp + A]. We then introduce the follow- 
ing fermion operators: 



Cfe,, 




folk e R 
for /c e L 
otherwise. 



(2.3) 
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kj 




-7C/a \ 










FIG. 1: Single-particle energy band. The annihilation operator of 
an electron near the Fermi points with momentum k £ [— A;f — 
A, — /cf+A] (fc e [/cf — A, fcF+A]) is denoted ak,-,tT (afc,+,o-), and 
that of an electron far away from the Fermi points is denoted bk^^r- 



The index \\ (_L) of the coupling constants denotes the scatter- 
ing of electrons with same (opposite) spins. 



B. Vertex corrections 



In the conventional weak-coupling approach to the ID 
EHMr^ii^ one estimates the coupling constants in Eq. i2.4\ 
only up to the lowest order in U and V: 



Electrons near the Fermi points are shuffled by the two- 
particle scattering: Hint ~ U J2j "i,T ^j-l + ^ J2j "-j+i- 
Following the standard g-ology approach, we will focus 
on the scattering processes between electrons near the Fermi 
points, i.e., the scattering processes which involve ak,±.a- only. 
The Hamiltonian for such interaction processes is 



Hi 



int 



.91 II 

ki ,p,(T 

2L ^ 



t 



ki .,p,a 



.9211 

2L 



ki,p,<j ^fe2,-p,o- ^k3,-p.a ^k4.,p,<j'- 



+ 



■'^ll,p,cr ^fc2:P,o- "/C3,-P,cr "fc4,-p,a-- 

ki^p,(7 

ff2-L . t t 

2^ '^ki,p,cr ^k2,p,<y '^k3,-p,a '^k4,,-p,a- 

ki ,p,a 

^'ki,p,<j '^fc2,-p,<T ^\3.,p,a '^k4,,-p,<j'- 



M .„t 

ki ,p,cr 

/ . •^ki,p,(T ^k2,-p,cr ^k3,p,a ^k^, 
ki ■,p,(y 

54|| t t 

2^ -al aj.„ „ „ aL a 

ki,p,a 



2L 



^p,a ■ 



54_L \ - t 

2L 2Z 

ki ,p,(T 



ki,p,tT k2,p,<T k3,p,a ki,p,a ■ 
t 



ki,p,<y ^k2,p,(7 '^k3,p,a ^ki, 



(2.4) 



where ct = i (t) for a = | (4), i is the length of the system, 
and : : denotes normal ordering. The summation over the mo- 
mentum ki is taken under the condition of the total momentum 
being conserved (equal to ±27r/a for the umklapp scattering). 
The index p = +/— denotes the right-Zleft-moving electron. 
The coupling constants gi|| and gi±^ {g^ and .gsx) denote the 
matrix elements of the backward (umklapp) scattering, while 
(/2|| and g2± {g^w and .g4j_) denote the matrix element of the 
forward scattering with the different (same) branch p — ±. 



911- 


= 93- 




= 9A- 


9i\\ 


= .931 


52 II 


= 94\ 



= {U~2V)a, 
= {U + 2V)a, 
= -2Va, 
= +2Va. 



(2.5a) 
(2.5b) 
(2.5c) 
(2.5d) 



In analyzing the low-energy physics of Eq. ( 12. 4> . one then 
employs the standard g-ology approach,^"* i.e., the perturba- 
tive RG method, and obtains flow equations for the marginal 
terms in Eq. ( 12. 4> . From this RG analysis'^-^^ one finds that 
the g3± term generates a gap in the charge excitation spec- 
trum if \g3±\ > -{g2\\ + g2± - 9i\\) and 53^ ^ 0, whereas 
the gi± term yields a gap in the spin excitation spectrum if 
\gi±\ > -(.9211 - .92± - .9i||) and gi± ^ 0. Hence, with 
the lowest-order coupling constants Eq. ( 12. 5> , one would con- 
clude that the charge (spin) excitations become massless at 
U - 2V = (U - 2V > 0). This would mean that, as 
U increases, both the charge and spin sectors become criti- 
cal simultaneously atU — 2V, where a direct and continuous 
CDW-SDW transition takes place. This analysis is found to be 
insufficient from the following argument. The (accidental) si- 
multaneous vanishing of g^^ and gi± results from the lowest- 
order estimate in U and V and there is no symmetry princi- 
ple that enforces gi± and g^^ to vanish simultaneously. It is 
possible that the higher-order corrections to g lift the degen- 
eracy of zeros and change the topology of the phase diagram. 
Therefore, in order to analyze the phase diagram at [/ w 2V, 
we need to go beyond the lowest-order calculation of the cou- 
pling constants in the g-ology. In this section, we compute 
the vertex corrections due to virtual processes involving high- 
energy states^ by integrating out bk,a- This procedure allows 
us to obtain the effective coupling constants g's that include 
higher-order corrections. 

The second-order vertex diagrams for the coupling con- 
stants are shown in Fig. |2] The solid lines denote the low- 
energy states ak,±,a, while the dashed lines denote high- 
energy states 6fc The nonzero contributions from the 
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A only weakly, and we can set A = 7r/4 in the following anal- 
ysis as we are interested in the qualitative feature of the phase 
diagram (different choices will only lead to small quantitative 
changes in phase boundaries). Incidentally, the logarithmic 
divergence of Ci (A) in the hmit A — > leads to the familiar 
one-loop RG equations. 



C. Bosonization 



FIG. 2: Vertex diagrams with second-order corrections [(a)-(e)]. 
Solid lines denote electron states in the momentum space k £ R 
or k £ L, while the dashed lines denote electron states in the other 
momentum space. 



Having integrated out the high-energy virtual scattering 
processes, we now focus on the low-energy states and lin- 
earize the dispersion of ak.±.cr around the Fermi points. The 
kinetic-energy term with the linearized dispersion is given by 



second-order virtual processes (a)-(e) are 



^921 

^9n 



^9% 



Ms 



u 



int nt 



nt 



nt 



D2a, 



(2.6a) 
(2.6b) 
(2.6c) 
(2.6d) 



T^T- Dia-—D2a, (2.6e) 

Ant nt 



' nt 



(2.6f) 



([/- 2^)2+41^2 y2 

^ -i-- D^a + —D2a, (2.6g) 

47rt nt 



where 



Di(A) ^ 
D2{A) = 



fix l2-aK 



-7r/2+aA 



dk 
cos k ' 



dk ■ 



7r/2+aA 



COS k 



(2.7a) 
(2.7b) 



By introducing Ci(A) = 2 ln[cot(aA/2)] and C2(A) = 
2 cos aA, Di (A) and D2 (A) are rewritten as Di (A) = Ci (A) 
and Da (A) = Ci(A) - C2(A). In terms of Ci and C2, the 
couphng constants with second-order corrections are given by 



gi^^iU-2V)a 
Ci 



^V^a, (2.8a) 
nt 



9i\\ 



2Va- —(U-2Vya 

Ant ^ ' nt 



(2.8b) 



53± = iU- 2V)a 
Ci 



G 



^V'a, (2.8c) 
nt 



g3^^=-2Va~^(U-2Vfa+^V\, (2.8d) 
Ant nt 

and g2\\ — +2Va, g2± = {U + 2V)a, 54j| = +2Va, and 
g4± = {U + 2V)a. Except when aA ^ 1, the C^s depend on 



VF{-k-kF)al^_^^ 



(2.9) 



keL,a 



where vp — 2ta is the Fermi velocity. The field operators of 
the right- and left-moving electrons are given by 



'^-,<y{x) 



Vl 



ikx 



keR 



(2.10a) 
(2.10b) 



keL 



We apply the Abelian bosonization method and rewrite the 
kinetic-energy term Hq — J dx Hq in terms of bosonic phase 
fields as (see AppendixlAl 



Ho 



VF 

An 



{2nne 



VF 

An 



(2.11) 



where (0) is the bosonic field whose spatial derivative is pro- 
portional to the charge (spin) density, [0{x), (f){y)] = 0. The 
operators Ilg and II^ are canonically conjugate variables to 9 
and (p, respectively, and satisfy the conventional commutation 
relations, [e{x),ng{x')] = [(j){x),n^{x')] = iS{x-x'). We 
also introduce chiral bosonic fields 



(j)±{x) 



(x) =F 2n 
{x) T 27r 



dx' ng{x') 
dx'U^{x') 



(2.12) 
(2.13) 



One can easily verify that these chiral fields sat- 
isfy the commutation relations [d±{x),0±{x')] = 
[4>±{x),(j)±{x')] = ±i(7r/2) sgn(a; — x') and 
[e+{x),e-{x')] = [(j)+{x),(j)-{x')] = in/2. In terms 
of these fields, the kinetic-energy density reads 



^0 = 5 E [(^^^p 



(2.14) 
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To express the electron field operators ^|Jp,a■ with the 
bosonic phase fields, we introduce a new set of chiral bosonic 
fields 

V'p.T = + 4>pi VpA ^ Sp- 4>p^ (2.15) 
which obey the commutation relations 

[^±,crix),ip±^a'ix')] = ±iTTSgn{x - x')S„,a', (2.16a) 
[ip+^crix),ip-,a'ix')] = iirSca'- (2.16b) 

In terms of the phase fields (pp^a, the electron field operators 
can be written as 

i'p,a{x) = —j== exp [ipkpx + ipifp.aix)] , (2.17) 
V27ra 

where the Klein factor rja satisfies the anticommutation re- 
lation {rjajVa-'} — 25cr^o-'. One can verify that the operator 
defined in Eq. ( I2.17> satisfies the same anticommutation rela- 
tion as the fermion field operator It follows from Eq. ( I2.17> 
that the order parameters in Eq. j2.2> are rewritten as 

OsDw{x) cx cos6'(a;) sm(f>{x), (2.18a) 

CcDw(a^) oc sin6'(a;) cos(f>{x), (2.18b) 

CBCDw(a;) oc cos6'(a;) cos(/)(x), (2.18c) 

Cbsdw(2:) oc sin6'(a;) sin4>{x). (2.18d) 

The interaction part of the Hamiltonian Hint, Eq. i2A\ . can 
be also expressed in terms of the boson fields 9± and (l)± . It 
has been suggested that, besides the marginal operators, op- 
erators with higher scaling dimensions can play an impor- 
tant role in the first-order CDW-SDW transition^^-^"* which 
is known to occur in the strong-coupling region of the ID 
EHM. - T-^Pi^Pt^i We thus include all the terms of scaling dimen- 
sion 4 [= 2 (charge sector) + 2 (spin sector)]. We also note 
that there are some complications and subtleties in bosonizing 
the off-site interaction term, i.e., the nearest-neighbor inter- 
action term V (see Appendix IaI for detail). We obtain the 
bosonized Hamiltonian density, 

^ = ^ E [Md.Opf + v.idxcl^pf] 
P=+. - 

--^cos2^^ + -^cos2^ 
2n''a^ 27r"^a^ 

— ^^'l'^ „ cos2^^ cos 20 
27r^a^ 

- i!^{d^e+){d^0_) cos2</) 

+ ^ {d^(l)+) (d,^^) cos 20 
27r^ 

+ {9x0+) {dj^) {dx4>-) ■ (2.19) 

The renormalized velocities are Vp ~ 2ta + (54 1| + g4± — 
(7i||)/27r and v„ = 2ta + {g^ — g4± — gi||)/27r. The marginal 
terms with the couplings gp and gc {ga and gs) determine 



low-energy properties of the charge (spin) modes,^- where 

9p = g2± +92\\ - 5l||' 9c = 93±, 9a- = 92± - 92\\ 

gs = 9i±. The g^s, gps, gca, and gp^y terms with scaling di- 
mension 4 couple the spin and charge degrees of freedom. The 
gcs coupling comes from the umklapp scattering g^ . The gps 
igpa) coupling is generated from the backward scattering of 
antiparallel- (parallel-) spin electrons while the gca coupling 
is generated from the umklapp scattering of electrons with an- 
tiparallel spins (see Appendix IaIi. These coupling constants 
are given by g^s = gps = 9ca = 9pa = -2Va to lowest or- 
der in V . Cannon and Fradkin examined the effect of the g^s 
term and argued that it plays a crucial role in the first-order 
CDW-SDW transition. Voit included the gps and gca terms, 
as well as the gcs term, in the perturbative RG analysis of the 
coupling constants, but did not consider the gpc term.^"^ Here 
we note that it is important to keep the gp^ term as well, since 
the global SU(2) symmetry in the spin sector is guaranteed 
only when g^ = gs, gcs = gca, and gps = g pa- 



ly. Renormalization-group equations 

We perform a perturbative RG calculation to examine the 
low-energy properties of the ID EHM in the weak-coupling 
regime, taking into account quantum fluctuations of the phase 
fields. The operator product expansion (OPE) technique al- 
lows us to systematically handle the higher-order terms in 
the bosonized Hamiltonian ( I2.19> . The one-loop RG equa- 
tions that describe changes in the coupling constants during 
the scaling of the short-distance cutoff (a ae'*') are given 
by (see AppendixfSlfor their derivation) 





+ 2 G^ + G^^, + Gs Gps , 


(2.20) 


jr- - 


+ 2 Gp Gc — Gs Gcs ^ Gcs Gps, 


(2.21) 


^G - 

jr- - 


/^2 /^2 


(2.22) 


'^G 

jr^' - 


~ 2 Gcs + 2 G p Gcs ^ 4 Gs Gcs 




— 2 Gc Gs 2 Gc Gps 4 Gcs Gps , 


(2.23) 


'^G - 


^ 2 Gps + 2 Gp Gs 






-^GcGcs-^Gls-AGs Gps, 


(2.24) 



where G^ are dimensionless coupling constants with the ini- 
tial values Gi/(0) — gi,/{4:Trta). The number of the inde- 
pendent coupling constants is five, since the SU(2) spin sym- 
metry guarantees the relations G^ = Gs, Gca = Gcs, and 
Gpa = Gps to hold in the scaling procedure. From these scal- 
ing equations, one finds that the Gp, Gc, and Gs terms are 
marginal (the scaling dimension=2),^''-^^ while the Gcs and 
Gps terms are irrelevant operators of the scaling dimension 4. 
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III. PHASE DIAGRAM OF THE HALF-FILLED 
EXTENDED HUBBARD MODEL 

A. Bond-charge-density-wave state 

In this section, we show that the BCDW phase exists in be- 
tween the CDW and SDW phases in the weak-coupling region 
of the IDEHM. 

First we focus on the weak-coupUng limit U,V <^t, where 
we can neglect the irrelevant terms of scaling dimension 4 and 
restrict ourselves to the marginal terms cx gp, g^, gc, and gg. 
Effects of the irrelevant terms are discussed later in this sec- 
tion. Within this approximation, the Hamiltonian reduces to 
two decoupled sine-Gordon models, and we can analyze the 
properties of the spin and charge modes, separately. The one- 
loop RG equations for these coupling constants are given by 



with Gc 



Gn 



0: 



2Gl{l), 

2Gp{l)Gc(l), 

-2Gl{l). 



(3.1) 
(3.2) 
(3.3) 



The spin excitations are controlled by the Gg coupling, 
which is marginally relevant (marginally irrelevant) when 
Gs < {Gs > 0). If .gs < 0, then |Gs(OI increases with 
increasing I. In this case the phase field (j) is locked at (j) = 
mod TT to gain the energy [see Eq. J2.19H . and consequently 
the spin excitations have a gap. On the other hand, if gs > 0, 
then |Gs(OI decreases to zero as / increases, and the ip field 
becomes a free field; the spin sector has massless excitations. 
The approach of Gg to zero is very slow (^1/ 0, and the 
field has a strong tendency to be near = 7r/2 mod tt. Al- 
though it eventually fails to lock the phase 0, the marginally 
irrelevant coupling still has an impact on low-energy prop- 
erties by giving rise to logarithmic corrections to correlation 
functions.^** 

The charge sector is governed by the two couplings Gc 
and Gp, whose RG flow diagram is of the Kosterlitz-Thouless 
type. Since gp = {U + QV)a > 0, Gc is a relevant coupling 
and always flows to strong-coupling regime, unless gc = 0. 
This means that Gc(/) has two strong-coupling fixed points, 
Gc(/) — > oo and Gc{l) —oo, depending on its initial value 
f/c > and gc < 0. As seen from Eq. (I2.19> . the relevant gc 
with positive (negative) sign implies the phase locking of 6 at 
the position 9 — {tt/2) mod tt. 

From the above standard arguments, the ground state can be 
identified by simply looking at the initial value of the coupling 
constants gc and gg. The ground state is classified into four 
cases as summarized in Table |I] and the positions of locked 
phases {9, </)) for respective cases are shown in Fig.|3] 

(i) gs < and gc < 0'. The phase fields are locked at 
{9, (j)) = ((7r/2) + 7r/i, 7r/2), where Ii and I2 are integers. In 
this case, among the order parameters in Eqs. (I2.18> . only the 
CDW order parameter has a finite expectation value, and the 
ground state is found to be the CDW state. Both charge and 
spin excitations are gapped. 



TABLE I: Possible ground-state phases and positions of (quasi) 
locked phase fields determined only from the marginal terms in Eq. 
<2T9l . 



Phase 




{9c,9s) 


SDW 


i0,±7v/2),(n,±n/2) 


(+,+) 


CDW 


(±V2,0),(±7r/2,^) 




BCDW 


(0,0),(7r,7r),(0,7r),(7r,0) 


(+,-) 


BSDW 


i±7v/2, ±n/2) 


(-,+) 



(ii) gs < and gc > 0: The phase fields are locked at 
{9,(f)) = (7r/i,7r/2)- The nonvanishing order parameter is 
then Obcdw. and the ground state is the BCDW state. Both 
charge and spin excitations are gapped. 

(iii) gs > and gc < 0'. The field 9 is locked at 6* = 
{tt/2) + nil, and the field cj) tends to be around (f> — (7r/2) + 
7r/2 although it is not locked in the low-energy limit. In this 
case the dominant correlation is that of the BSDW state. The 
charge excitations are gapped whereas the spin excitations are 
gapless. 

(iv) gs > and gc > 0: The field 9 is locked at 9 = tt/i, 
whereas the field cj) tends to be near (j> = {tt/2) + ttI2- The 
dominant correlation is the SDW order The charge excita- 
tions are gapped while the spin excitations are gapless. 

Combining the results of Table|I]and the coupling constants 
Eqs. ( 12. Sai l and ( 12. Sel l, we obtain the ground-state phase dia- 
gram of the ID EHM in the weak-coupling limit. For U larger 
than 2V such that gc > and gs > 0, we have the SDW phase, 
while for U sufficiently smaller than 2V (gc < and < 0) 
we have the CDW phase. At U = 2V, we see from Eqs. 
( 12. Sat and ( 12. Set that gs{~ gi±) < and gc{— 53_l) > due 
to the G2 term. This implies that a new phase different from 
the CDW and SDW states appears for U w 2V. From Table|I| 
we identify the new phase with the BCDW phase. Within the 
approximation we employ here, the phase boundary between 
the BCDW phase and the CDW (SDW) phase is located at 




• 


SDW 





CDW 


A 


BCDW 


▲ 


BSDW 



FIG. 3: Positions of locked phase fields 9 and in the SDW, CDW, 
BCDW, and BSDW states. 



gc = (gs = 0). In this phase diagram, the charge excitations 
are gapful except on the CDW-BCDW transition line, while 
the spin excitations are gapless in the SDW phase and on the 
SDW-BCDW transition line. From Eqs. ( l3.H -( l33t . we can 
estimate the charge gap Ac and the spin gap As as 

Ac^M^j , As^texpi—-] (3.4) 

for \gc\ ^ gp ^ ta and < —gs ^ ta, respectively. 

Next we examine effects of the parallel-spin umklapp scat- 
tering gcs on the BCDW state. We consider the situation very 
close to the CDW-BCDW transition by assuming gc ~ and 
gs < 0, i.e., U~2V ^ -CzFVTrt + OiV^/t"^). In this case 
the spin gap is formed first as the energy scale is lowered. For 
energies below the spin gap, we can replace cos 20 with its 
average (cos 20) « (A^/t)^. This means that the coupling 
constant gc is modified as 

9c 9c + gcs (cos 2(f)) . (3.5) 

Thus we find that the BCDW state, which is realized for 
g* > 0, becomes less favorable due to the gcs{< 0) term. 
We note, however, that the CDW-BCDW boundary does not 
move across the U = 2V line because j^cs (cos 20)1 « 
2Vaexp[—c{t/V)'^] is much smaller than the C2 term in Eq. 
(I2.8c> for V <^ t, where c is a positive constant. A similar 
argument applies to the region near the SDW-BCDW transi- 
tion. Suppose that U-2V ^ +C2V^/7rt + 0{V^/t^) where 
gs w and gc > 0. In this case, as the energy scale is low- 
ered, the charge gap opens first and the 6 field is pinned at 
6 = mod TT. Below the charge-gap energy scale, the field 
is subject to the pinning potential g* cos 20 with 

g*s = 9s- gcs{cos2e), (3.6) 

where (cos 26*) w (Ac/t)^'^"'^''). Thus the BCDW phase, 
which is now realized for g* < 0, also becomes less favorable 
by the ~gcs{cos 26) (> 0) term. Again the phase boundary is 
not moved beyond the U = 2V line since | gcs (cos 2^?) | « 
2Va{ c'V/t)''^/^ is much smaller than the C2 term in Eq. 
(I2.8a> . where c' is a constant of order 1 . Therefore we con- 
clude that the BCDW phase is robust against the gcs term in 
the weak-coupling limit. The analysis in this section estab- 
lishes the existence of the BCDW phase near U w 2V for 

B. First-order SDW-CDW transition 

In this section, we discuss how the BCDW phase becomes 
unstable at strong coupling and how the two continuous tran- 
sitions change into the first-order SDW-CDW transition. 

To our knowledge. Cannon and Fradkin were the first to 
argue that the g^ term (describing the umklapp scattering of 
parallel-spin electrons), which is conventionally ignored due 
to its large scaling dimension, can become relevant at large 
U and V and cause the first-order CDW-SDW transitioni2^ To 
get an insight into the effect of the gcs term in the relevant 
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FIG. 4: Phase diagram obtained by minimizing the potential V{0, 0) 
for gcs < 0. The double line denotes the first-order transition, while 
the single line denotes the second-order transition. Bicritical points 

are at {gc,gs) = {±\gcs\,T\gcs\)- 

case, we perform a semiclassical analysis: we neglect spatial 
variations of the fields in Eq. ( I2.19> and focus on the potential, 

V{6, 0) = —gc cos 26 + gs cos 20 — gcs cos 26 cos 20, (3.7) 

where gcs — g3\\ < 0. The order parameters of the 
SDW, CDW, BCDW, and BSDW states take maximum am- 
plitudes when the fields 6 and are pinned at {6, 0) = 
(7r/i,(7r/2) + TTh), ((7r/2) + nluTrh), {nlunh), and 
((7r/2) + nil, ('''/2) + 71/2), respectively, where Ii and I2 
are integers. The potential energy in these states is obtained 
by inserting these pinned fields into Eq. ( I3.7> . e.g, Vsdw = 
V{7rli, {tt/2) + Tr/a), yielding 

Vsdw = -gc - gs ~ \gcs\, (3.8a) 

VcDw = +gc + gs - \gcs\, (3.8b) 

Vbcdw = -gc + gs + \gcs\, (3.8c) 

Vbsdw = +gc- gs + \gcs\- (3.8d) 

We find that the gcs term stabilizes the SDW and CDW states 
while it works against the BCDW and BSDW states. Compar- 
ing these energies, we obtain the phase diagram in the gc-gs 
plane at a fixed gcs (Fig.l^. In the presence of the gcs term, 
the direct CDW-SDW transition line appears in this phase di- 
agram. 

We now discuss the nature of the phase transitions. The 
potential V{6, 0) on various transition lines is shown in Fig. 
|5] On the boundary between the SDW and BCDW phases, 
which is located at gs = —\gcs\ and gc > \gcs\, the potential 
takes the form V{6, 0) = —gc cos 26 + gs cos 20(1 — cos 26) 
[Fig-EIa)], which pins the 6 field at 6* = tt/i and leaves the 
field completely free. We thus find that the SDW-BCDW 
transition is continuous, i.e., the SDW and BCDW phases co- 
exist without potential barrier on the phase boundary. On the 
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FIG. 5: The potential V{e, (j)) on the SDW-BCDW (a), BCDW- 
CDW (b), and CDW-SDW (c) transition Unes. 



boundary between the BCDW and CDW phases, located at 
5c = \9cs\ and gg < ~\gcs\, the potential now takes the 
form V{9,(p) — ~gcCos26{l — cos 20) + gsCos2(/) [Fig. 
|5jb)]. The potential locks the (p field at = 71/2, where it 
has no effect on the 9 field. Thus, we find that the CDW- 
BCDW transition is also continuous. From similar consider- 
ations, we find that the SDW-BSDW and BSDW-CDW tran- 
sitions are continuous as well. In Fig. |3 the phase bound- 
aries of continuous transitions are shown by the solid lines. 
On the contrary, the phase boundary shown by the double 
line in Fig. |3 is of different nature from the others. The 
potential V{9, (f) on the double line is shown in Fig. |5jc). 




FIG. 6: The potential V{9,(f)) on the bicritical point ((/c,5s) = 
(Iffcs I, — |<?cs I). The potential minima are the lines 6 = tt/i and 
(j> — 7r/2. 

where the potential minima are given by the isolated points 
(9,(1)) = (tt/i, (7r/2) + nh) and {{n/2) + 71/1,71/2). These 
minima correspond to the SDW state and the CDW state, see 
Fig-El The point to note is that there is a finite potential barrier 
of height min ( | gcs | , 2 1 gc s | — 2 1 | ) between the corresponding 
minima for the SDW and CDW phases. Hence we conclude 
that the CDW-SDW transition is first order when gcs is rele- 
vant. 

From the above arguments, we find that strong umk- 
lapp scattering of the parallel-spin electrons destabilizes the 
BCDW and BSDW states and gives rise to bicritical points 
(.9c, <?s) — ±(5cs, ~.9cs) where the two continuous-transition 
lines merge into the CDW-SDW first-order transition line. Let 
us take a closer look at these bicritical points. Taking into 
account the fact that gc > and gs < for U w 2V in 
the original EHM, we will focus on the bicritical point at 
{gc,gs) = (Iffcsl, — IffcsD- The effective potential at the bi- 
critical point takes the form 

V{9,(j)) = -5(cos26' + cos2(/)-cos26'cos2(/i), (3.9) 

which is shown in Fig. |5] This potential has an interesting 
feature that its potential minima are not isolated points but 
the crossing lines 9 — nm or (j) — nn (m, n: integer). On 
these lines either 9 or (j) becomes a free field; the theory has 
more freedom than a single free bosonic field, but less than 
two free bosonic fields. We thus expect that the theory of the 
bicritical point should have a central charge larger than 1 but 
smaller than 2. Detailed analysis of the critical theory is left 
for a future study. We note that when gcs = the first-order 
CDW-SDW transition line collapses into a tetracritical point, 
{gc,9s) = (0,0), and the phase boundaries in Fig.l^reduce 
to the lines gc — and gs ^ where all the transitions are 
continuous. 

Fabrizio et alM and Bajnok et alw^ discussed effects of 
higher-frequency terms, such as sin 36* and cos 46*, which are 
generated through the renormalization-group transformation. 
From the semiclassical arguments, it can be seen that these 
terms can also change a second-order transition to a first-order 
transition. In fact, it was argued that these higher-frequency 
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FIG. 7: Phase diagram of the half-filled ID extended Hubbard 
model. The double line denotes the first-order transition, while the 
single lines denote the second-order transitions. The bicritical point 
is at ((7c, 14) ~ (5.0t, 2M). 



terms make the SDW-CDW transition first order in the strong- 
coupHng regime of the ID EHM.^^ However, we have shown 
that the SDW-CDW first-order transition can occur simply due 
to the gcs term which is the leading irrelevant term in this sys- 
tem. Since the higher-frequency terms are even less relevant 
than the g^s term, we expect that the g^s term should play a 
dominant role in the first-order transition in the ID EHM. 



C. Global ground-state phase diagram 

To obtain the global phase diagram of the ID EHM, we 
have numerically solved the scaling equations ( I2.20l l-( |2.24| |. 
We find out which phase is realized by looking at which one 
of the couplings Gc, Gg, and Gcs becomes relevant first, as 
we have discussed in Sees. IIII Al and 1111 B I First, if \Gc\ 
grows with increasing I and reaches, say, 1 first among the 
three couplings, then we stop the integration and compute 
G* = Gs — GcsSgn(Gc). Since the charge fluctuations are 
suppressed below this energy scale, we are left with Eq. j3.3> . 
where Gg is replaced by G* . We immediately see from TableU 
that a positive (negative) G* leads to the SDW (BCDW) state 
for Gc > and the BSDW (CDW) state for Gc < 0. Second, 
if \Gs \ becomes 1 first, or more precisely, if Gg reaches — 1 
first, then we are left with Eqs. ( 13. 1> and ( 13. 2K where Gp and 



Gc are replaced by G* ~ Gp 



Gps and G* 



Gc + Gc 



respectively. We see that a positive (negative) G* leads to 
the BCDW (CDW) state. Finafly, when |Gcs| reaches 1 first, 
we stop the calculation and compare Gc and Gg- Since both 
charge and spin fluctuations are already suppressed by the 
Gcs cos 29 cos 2(p potential, we can deduce the phase from the 
semiclassical argument. From Fig.0]we see that we have the 
SDW state for Gg > -Gc and the CDW state for Gg < -Gc- 
Here we note that in the SDW state the pinning potential to the 



(f) field is marginally irrelevant and thus the spin sector should 
become gapless. 

The phase diagram obtained in this manner is shown in 
Fig. The single lines denote continuous transitions, and 
the double line denotes the first-order transition. In the weak- 
coupling limit, the BCDW phase appears atU ^ 2V and the 
successive continuous transitions between the SDW, BCDW, 
and CDW states occur as V/U increases. When U and V 
increase along the line U w 2V, the BCDW phase first ex- 
pands and then shrinks up to the bicritical point {Uc,Vc) ~ 
{5.0t, 2.3t) where the two continuous-transition lines meet. 
Beyond this point the BCDW phase disappears and we have 
the direct first-order transition between the CDW and SDW 
phases. The phase diagram (Fig.0 is similar to the ones ob- 
tained by using more sophisticated numerical methodsJi^^ 
We note that the position of the first-order transition line in 
Fig. Q is not reliable quantitatively as we have used the per- 
turbative RG equations. The recent Monte Carlo calculation?^ 
gives the most reliable estimate for the position of the bicrit- 
ical point, (C/c, Vc) ~ ((4.7 ± 0.l)t, (2.51 ± OM)t), which 
agrees with our estimate in Fig. within 10%. The semi- 
quantitative agreement gives us confidence that our approach, 
semiclassical analysis of the low-energy effective Hamilto- 
nian derived with use of the perturbative RG, is reliable even 
in the strong-coupling regime near the multicritical point. 



IV. EFFECT OF STAGGERED SITE POTENTIAL 

In this section, we examine effects of alternating on-site 
modulation of the chemical potential, i.e., the staggered site 
potential, in the half-filled ID EHM. The Hamiltonian to be 
considered is given hy H' = H + Ha with H defined in Eq. 
dO and 



i/A = A^(-l) 



(4.1) 



The model is called the ionic Hubbard model if V = 0. 
When U =^ V = 0, the system is a trivial band in- 
sulator, since the A term induces a gap 2|A| at k ~ 
±7r/2 in the single-particle spectrum and the lower band is 
fully filled. For many years effects of the on-site repul- 
sive interaction U on the band insulator have been investi- 
gated intcn:i-?l y^^-^^-^'*-^^-^^-^^-^^-^^-'^°-"^'-'*^-'*^-'^'^-'^^-'^^-'^^-'^^-'^^-^°-^^ 
from both numerical and analytical approaches. Using the 
standard bosonization method, Fabrizio, Gogolin, and Ners- 
esyan recently argued that the ground state of the ionic Hub- 
bard model exhibits three phases as U increases: the band in- 
sulator, the SDI, and the Mott insulator. The order parameter 
of the SDI state is nothing but that of the BCDW state, and we 
can regard the two states as essentially identical. It was also 
argued that the quantum phase transition from the band insu- 
lator to the SDI state belongs to the Ising universality class 
whereas the other transition from the SDI state to the Mott in- 
sulator is of the Kosterlitz-Thouless type. Recent numerical 
studiespi^Si^^i^i^2i^2ii2iiiilS however, have reported controver- 
sial results on the existence of the SDI phase. Some claimed 
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to find two quantum phase transitions while others found ev- 
idences of only one phase transition. With this issue of the 
SDl phase in mind, in this section we investigate the phase di- 
agram of the ID extended Hubbard model with the staggered 
site potential and examine critical properties of the quantum 
phase transitions. 

We take into account the staggered site potential and the 
correlation effects on equal footing by treating them as weak 
perturbations. We use Eq. ilAli to rewrite Ha in the contin- 
uum limit as Ha — J dx Ha, whera^^i^ 



Ti-A 



3A 



2{TTay 



sin 6* cos( 



(4.2) 



with gA = 47rAa. Note that the CDW order parameter Ocdw 
is proportional to Ha, and gA can be regarded as an exter- 
nal force coupled to Ocdw- This has the consequence that 
Ocdw acquires a nonvanishing expectation value for any fi- 
nite U and V, as long as gA ^ 0. In this section we will de- 
note the insulating phase connected to the free-electron band 
insulator (U — V ^ and A ^ 0) by the Bl phase, rather 
than the CDW phase. 

The bosonized form of the Hamiltonian H' can be thought 
of as a generalization of the so-called double sine-Gordon 
(DSG) model as H' contains sine/cosine terms with differ- 
ent frequencies (sin 6* and cos 26*, cos0 and cos 20). The 
DSG theory itself has been investigated intensively^2iS2iffi and 
shown to have a critical point belonging to the Ising univer- 
sality class [c = i conformal field theory (CFT)]. To obtain a 
qualitative understanding of the critical properties in our sys- 
tem, we first perform a semiclassical analysis in a similar way 
to Sec. Ill B, before examining the global phase diagram of 
H' with use of the RG method. 



A. Semiclassical analysis 

In this section, we perform a semiclassical analysis to the 
Hamiltonian H' = H + Ha, where H and Ha are given by 
Eqs. ( I2.19> and i4.2l . respectively. We neglect spatial varia- 
tions of the field and focus on the locking potential: 



VA{0,(j)) = —gcCos29 + gsCos2(f> — 
— gA sin 9 cos 0. 



gcs cos 26 cos 2(j) 
(4.3) 



First, we examine the case gcs = 0, which corresponds to 
the situation where the gcs term becomes UTelevant in the RG 
scheme. The potential to be considered is 

^ VA{e,(i>)\g^^=o 

= —gc cos 29 + gs cos 2(f) — gA sin 9 cos (j). 

(4.4) 



Vl{9,< 



Due to its double-frequency structure, possible locations of 
the phase locking are different from the ones we found in 
Sec. Ill B. For example, when gc > (.gs > 0), the two 
kinds of potentials proportional to sin 9 and cos 29 (cos cj) 
and cos 2(f)) compete with each otheri^M£ The locking of the 



TABLE II: Possible ordered ground states and the position of (quasi- 
)locked phase fields determined from Eq. 14. 4> . 



Phase 



SDW 

BI (for gA > 0) 
BI (for gA < 0) 
BCDW 
BSDW 



(0,±7r/2),(7r,±7r/2) 
{+n/2,0),{-7T/2,n) 
{+7T/2,n),{-n/2,0) 
{+{7T/2)±al0),{-{n/2)±aln) 
{+n/2,0±a°),{-7T/2,±{n-al)) 



phases 9 and are determined from the saddle-point equa- 
tions: cos 9{Agc sin 9 — gA cos 0) = and sin (f){~4gs cos (f + 
gA sin 9) = 0. In order to simplify the notations, let us intro- 
duce 



,0 _ 



cos 



5A 



cos 



9A 

4(7s 



(4.5) 



< TT are 



where \gA/gc\ < 4, \gA/gs\ < 4, and < aO,aO 
assumed. The solutions of the saddle-point equations yield the 
following four states with distinct configurations of the locked 
phase fields 9 and (j) (modulo 27r): (i) the SDW state with 9 
and locked at {9, (j)) = (0, ±7r/2) or (tt, ±7r/2); (ii) the BI 
state with {9, (p) = (+7r/2, 0), (-7r/2, tt) if gA > and with 
(6», 0) = (+7r/2, tt), (-7r/2, 0) if < 0; (iii) the "BCDW" 
state where the BCDW order and the CDW order coexist and 
which is realized when (9, 0) — (7r/2 ± a^, 0) or {—tt/2 ± 
a°,n); (iv) the "BSDW" state where the BSDW and the CDW 
order coexist and which is realized when {9, </>) = {tt/2, ± 
a^) or (— 7r/2, ±(7r — a°)). Tablelllland Fig.|8]summarize the 
possible ordered ground states and corresponding positions of 
locked phase fields. The potential energies in these states are 
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FIG. 8: Positions of locked phase fields 9 and (j) in the four states 
when gA > 0- 
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the form 



[sin — sgn((7A) cos 0] } 



BCDW 



FIG. 9: Phase diagram obtained by minimizing the potential energy 
V2.{0, <f)) [Eq. <4.4H . The phase boundaries between the SDW state 
and the BCDW state, and between the SDW state and the BSDW 
state are given by the curve Qs = g\l{^&gc) with Qc > 0. The phase 
boundaries between the BI state and the BCDW state, and between 
the BI state and the BSDW state are given by the lines gc = j \gA \ 
with gs < I jgAl and = 3 Iff a] with gc < j\gA\, respectively. 
All the phase transitions in this figure are continuous. The tetracriti- 
cal point is located at (ffcffs) = {j |5a|, j |ffA|). 



given by 



1/0 
•^SDW 

I/O 

"bcdw 



I/O 

•'BSDW 



= -.9c 



+9c + gs - \gA\ 
9a 

-gc + cjs - -J—, 
8.gc 

gl 

+gc- gs - -3— • 



(4.6a) 
(4.6b) 

(4.6c) 
(4.6d) 



In deriving Eqs. ( I4.6ct and ( I4.6dt . we have assumed 
IffA/ffcl < 4 and |(?a/3s| < 4, respectively. The CDW state 
is stabilized strongly by the term whereas the BCDW state 
and the BSDW state are also stabilized by the second-order 
contribution of ^a- By comparing these energies, we arrive 
at the phase diagram shown in Fig. |9] As we go across the 
boundary (gc — \ gA) from the BI state to the BCDW state, 
we find that each potential minimum splits into two minima, 
e.g., (0, 0) = (7r/2, 0) ((vr/2) ± a°g, O), and that the po- 
tential for the 9 phase field takes a double-well structure in 
the BCDW state. Similarly, as we go from the BI state to the 
BSDW state, each potential minimum splits into two minima, 
e.g., (0,0) = (7r/2,0) (7r/2,±a^), and now the poten- 
tial for the (f) phase field has a double-well structure in the 
BSDW state. As long as gcs = 0, any quantum phase transi- 
tion is continuous since a potential barrier between two poten- 
tial minima corresponding to two different states vanishes at 
the transition. The phase diagram (Fig.|9j indicates that a di- 
rect transition from the SDW state to the BI state takes place 
only when the parameters gc and gs are on the multicritical 
point (5c, 5s) = (j\gAl\\gA\), where the potential takes 



= ±(|7r 



if 



= ^|5a|{-1 
and is minimized at = ±[(7r/2) — 6*] and 
if 5A > 0, or at = ±[{tt/2) + 9] and : 
5A < 0. 

Let us take a closer look at low-energy excitations in the BI 
state and the BCDW state. The massive sine-Gordon model 
has topological excitations, solitons, and antisolitons. They 
are characterized by the topological charges Q and Sz for the 
charge and the spin sectors. 



= — I dxdrr.9, 

n 



2n 



dx 



(4.7) 



In the noninteracting case (U = V = 0) with a finite A, 
the lowest-energy excitation is a soliton of 9 and (p con- 
necting two neighboring minima of the — 5a sin 9 cos 0, e.g., 
{0,^)U^-^ = {-7T/2,n) and (0,0)U^oo = in/2,0). 
Such an excitation carries the charge Q ±1 and the spin 
Sz = which is nothing but a single-electron excita- 

tion in the band insulator. It has been pointed ouiSii^ii 
that in the SDl phase (i.e., in the BCDW phase), the topo- 
logical charge Q of the lowest-energy excitation becomes 
fractional, Q — ±2ag/TT, reflecting the local double-well 
structure of the potential near the potential minima, e.g., at 
[9, (j)) = [■K/2±a% 0). This is a unique feature of the BCDW 
phase and is contrasted from the integer charge Q — ±1 of 
the lowest-energy excitation in the pure BCDW phase where 
the phase fields are locked at [9, (f) = (0, 0). Accordingly, the 
phase transition between the BCDW state and the BI state be- 
longs to a different universality class from the one between the 
pure BCDW state and the CDW state discussed in Sec. 1111 BI 
In the former case, a small potential barrier in a double-well 
potential in the BCDW state vanishes at the critical point and 
the effective theory for the low-energy excitations is the "(ys"*" 
theory known to describe the Ising phase transition, rather 
than the Gaussian theory that governs the transition between 
the BCDW and CDW phases. 

One might expect that a similar semiclassical analysis can 
be applied to the spin field 0. Within the semiclassical ap- 
proach the topological charge Sz in the BSDW phase of Fig. 
Intakes a fractional value, ±a^^/{2'K). However, since the 
Hamiltonian has the global SU(2) spin-rotation symmetry, the 
SDW state and the BSDW state cannot have a true long-range 
order This implies that the phase field (j) cannot be localized 
except in spin-gap phases where (j) is locked at (</)) = mod 
TT. The global SU(2) symmetry thus prohibits the Ising crit- 
icality in the spin sector. In fact, the BSDW phase in Fig.|9l 
turns out to be just the BI phase. 

Let us now consider the situation in which gcs 7^ 0. In this 
case, the phase fields 9 and are locked in a similar way to 
the case gcs = 0, but and are modified into — > ae 
as, where 



and 



ae 



cos 



5A 



4(5c- \9cs\) 
5A 



4(.9s - \gcs\) 



(4.8a) 
(4.8b) 
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SDW 




FIG. 10: Phase diagram obtained by minimizing the potential en- 
ergy Va(S, <;A) [Eq. <4.3H for gcs < 0. The phase boundaries 
are given by gs = —gc + ^ |<?a| between the SDW and the 
BI states, g^ = \gcs\ + 3 IsaI between the BI and the BCDW 
states, gs — \gcs \ + j |3a| between the BI and the BSDW states, 

gs ~ —\gcs\ + 3a/[16(3c — Iffcal)] between the SDW and the 

2 



BCDW states, and gc = — Iflcs 
the SDW and the BSDW states. 



+ 5'a/[16(5s - \9cs\)] between 
Multicritical points are located 



at {gc,gs) = {+\gcs\ + j\gA\,--\gcs\ + j\gA\) and {-\gcs\ + 
J I 5a I, +|<?ca| + J IsaI). The single lines denote second-order tran- 
sitions, while the double lines denote first-order transitions. 



Here we have assumed Ig^/ {gc — \9cs\)\ < 4 and I^a/ {gs — 
\gcs\)\ < 4. The potential energies in the four states become 



^SDW — 

Vbi = 

^BCDW = 

^SDW = 



-9c - gs 
-gc + gs 



\gcs\, 
\gcs\ - 



|5A| 



-gc + ffs + \gcs\ 



9l 



-gc 



Kgc - l.gcs|)' 
2 

9a 

K9s - \gcs\)' 



(4.9a) 
(4.9b) 

(4.9c) 
(4.9d) 



By comparing these energies we obtain the phase diagram 
(Fig. I10> . In the limit (/a — > this phase diagram reduces 
to Fig. 0] One can easily find that the gcs term favors the 
SDW state and the BI state over the BCDW state and the 
BSDW state. The direct SDW-Bl transition line acquires a 
finite length in the phase diagram, like in Fig. |3 The analy- 
sis of critical properties of each quantum phase transition is 
more complicated than that in Sec. Ill due to the presence 
of two kinds of charge-spin coupled terms, the g^ and gcs 
terms. Along the phase boundary between the SDW state 
and the BI state, the potential energy is minimized at discrete 
points, {e,(j)) = (-7r/2,7r), (0,±7r/2), (7r/2,0), (7r,±7r/2) 
for 5A > 0, or at (6*, = (-7r/2,0), (0,±7r/2), (7r/2,7r), 
(tt, ±7r/2) for g^ < 0. These points correspond either to 
the SDW state or to the BI state (see Table HI}. Since any 
path connecting these potential minima has to go over a po- 
tential barrier, the direct SDW-Bl transition is first order In 



addition, both the transition between the SDW state and the 
BCDW state and that between the SDW state and the BSDW 
state become first order when gcs 7^ 0. On the phase bound- 
ary between the SDW state and the BCDW state, the poten- 
tial has isolated minima at {0,(1)) ~ (0,±7r/2), (tt, ±7r/2), 
(— 7r/2 ± ae,Tr), and (+7r/2 ± a0,O). The pinning of the 
phase fields at these minima corresponds either to the SDW 
state or to the BCDW state (see Fig. |8|i. On the multicriti- 
cal points at {gc,gs) = {+\gcs\-^j\gA\,-\gcs\-\-j\gA\) and 



{~\gcs\ 



3I.9AI 



I (7A I ), the potential takes the form 



Vl\e,<j)) = -|5cs|(cos26' + cos20-cos26' cos2(/)) 

+ ^ IsaI |-1 + [sin 6* - sgn(5A)cos0]^| , 

(4.10a) 

Vl^{e,(j)) = +|.gcs|(cos26' + cos2(/) + cos26' cos2(/)) 

+ ^ |.gA| |-1 + [sin 6* - sgn(5A)cos0]^| , 

(4.10b) 



respectively. The potential minima of V^^^(6', </>) and 
Vl^{e,(t)) are located at (61,0) = (-Tr/2,7r), (0,±7r/2), 
(7r/2,0), and (7r,±7r/2) for gA > and at (6*, 0) = 
(-7r/2, 0), (0, ±7r/2), {n/2, tt), and {tt, ±tt/2) for gA < 0. 



Finally, we note that even in the SDW state (the Mott insu- 
lator) the CDW order parameter has a nonvanishing expecta- 
tion value. This is because the alternating site potential Ha 
has the same form as the CDW order parameter Ocdw oc 
sin 6 cos (f). Even though the semiclassical analysis indicates 
that the phase fields are pinned, say, at {6, (j>) = (0, ±7r/2), 
quantum fluctuations of the fields around the pinning position 
lead to a nonvanishing (Ocdw)- This can be easily seen in 
the limit of small A, where 



{O 



CDW/ 



cx Tr 



exp 



-J dx{n + nA) 



sm a cos ( 



oc gATr 



exp 



dxH. SIT? 9 cos^ 



^0. 
(4.11) 



B. Renormalization-group analysis 



We perform RG analysis to take into account quantum fluc- 
tuations that are ignored in the semiclassical analysis. As in 
Sec. Ill, we obtain the RG equations using the OPE method 
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(see AppendixlEl: 

— Ga = + Ga + 2 ~ '-^^ '-^c 

3 3 3 

— Gp = + — Ga + 2 G^ + G^^ + G^ Gps, (4.13) 
^Gc ~ — — Ga + 2 Gp Gc — G<; Gcs ~ Gcs Gps, (4.14) 

— Gs= — — G\ — 201 ~ GcGcs ~ G^csi (4.15) 
al 4 

Ti^cs — ~ -7 Ga ^ 2 Gcs + 2 Gp Gcs — 4 Gs Gcs 

a( 4 

— 2GcGs ~ 2GcGps ~ 4:GcsGps, (4.16) 

^Gps = — — G^ — 2 Gps + 2 Gp Gs 

-4GcGcs-4GL-4GsGps. (4.17) 

The initial value of Ga(0 is given by Ga(0) = A/t, while 
those of the other coupling constants are given by G,y(0) — 
g^/{A-Kta). Since the RG equations are invariant under the 
sign change of Ga (Ga ^Ga), we can assume Ga(0) > 

without losing generality in the following arguments. 

We determine the ground-state phase diagram in a similar 
way as in Sec.|III] That is, we integrate the scaling equations 
(I4.12> -( |4.17> numerically and find which one of the couplings 
[Ga(0' Gc{l), Gs{l), and Gcs(/)] becomes most relevant. By 
doing so, we have encountered the following four cases. 

(i) The case where Gc{l) grows fastest and becomes 1 at 

1 — lp+. Below this energy scale (i.e., I > lp+), the charge 
fluctuations are suppressed and the phase field 6 is locked at 
6* = or TT. For the discussion of the ground-state properties 
we may first neglect the gA term since (sin 9) cos = 0. The 
Hamiltonian density Ti' then reduces to 

= g$](a,0p)2-!^G:(a.0+)(a,</>_) 
p 

+ ^G:cos2(/), (4.18) 
7ra 

where G* — Gs{lp+) ~ Gcs{lp+)- We immediately see that, 
if G* > 0, the spin excitations are gapless and the ground 
state is the SDW state. On the other hand, if G* < 0, then 
the operators proportional to G* are relevant [Gl{l) — > —00 
under scaling] and the phase fields are locked as {0, cf)) — 
(0, 0), (0, tt), (tt, 0), (tt, tt), which coiTesponds to the BCDW 
state with ag 7r/2 (i.e., ^a ^ 0), see Table HII This would 
become the BCDW state with ag < 7r/2 in a more realistic 
treatment where the term is not simply ignored. 

(ii) The case where |Gc(OI grows most rapidly and 
Gc{l) — 1 at Z = /p_. The phase field 9 is then locked 
at 9 — ±7r/2 for I > Ip^. Below this energy scale one 
can replace the sin 6* potential by its averaged value, i.e., 
smO (sinf?) = ±1. The effective Hamiltonian at I = Ip^ 



is given by 

ZTT ^ — ' TT 

P 

T ^G^ cos(/.+ -^G: COS20, (4.19) 
7ra^ 7ra^ 

where G*^ = G^Q-p-) and G* = Gs{lp-) + Gcs{lp-), and 
the sign — /+ of the Ga term corresponds to the position of 
the phase locking 9 = +(7r/2)/ - (7r/2). When G* > 0, 
the two G* terms are marginally irrelevant, and the only rel- 
evant operator is =F cos (p. Then the phase field (p is locked at 
(/) = or TT, depending on the position of the charge phase 
locking 9 = +(7r/2) or — (7r/2). On the other hand, when 
G* < 0, both GJ^ and G* terms become relevant. However, 
these terms do not compete with each other The only effect of 
the G^ term is to lift the degeneracy between the neighboring 
minima of — cos 20, and hence the position of the phase lock- 
ing is the same as in the case G* > 0. Therefore, regardless of 
the sign of G*, the resultant phase is found to be the Bl state 
with the phase locking at {9, (p) = {tt/2, 0) or (— 7r/2, tt). 

(iii) The case where either |Gcs(OI or |Ga(OI is most rel- 
evant. Then both charge and spin fluctuations are suppressed, 
and the classical treatment is sufficient at lower energy scale. 
In this case, we find to which phase the ground state belongs 
by substituting the parameters Gc{l) and Gs{l) into gc and gs 
inFig.[lO| 

(iv) The case where Gs{l) is most relevant and becomes 
— 1 at Z = Icr- Below this energy scale the spin fluctuations 
are suppressed and the phase field (f> is locked as — > or tt 
for I > Ic- The effective Hamiltonian of the remaining charge 
sector is 

= ^T.(W + ^G;id^9+){d,9^) 
p 

T^Glsm9~^G*c cos 261, (4.20) 

where G; = Gp{l„) - Gps(Z.), G\ = Ga(/.), and G^ - 
Gc{la) + Gcs{1<t)- The sign — /+ of the G\ term corresponds 
to the position of the phase locking </> = O/tt. In this Hamilto- 
nian, both of the nonlinear terms, sin 9 and cos 29, are relevant 
operators. If G* < 0, then the situation is the same as the case 
(ii): the G\ and G* terms do not compete with each other and 
the possible phase locking pattern is 6* = +7r/2 (— 7r/2) for 
(f) — Q (tt), where the ground state is the Bl state. If G* > 0, 
these two terms compete with each other, since the — (+) sin 9 
potential tends to lock the phase field 9 ai9 — +7r/2 (— 7r/2), 
while the cos 29 potential tends to lock it at = or tt. In this 
case, possible ground states are the Bl state and the BCDW 
state, and the quantum phase transition between them is of 
the Ising transition type with the central charge c = 1/2, as 
discussed in the preceding section. However, it is hard to esti- 
mate quantitatively the critical value of the coupling constants 
at the quantum phase transition. One way to estimate it is to 
find a critical point separating the basins of attraction to the 
two strong-coupling fixed points, (G^,G*) (+00,— 00) 
and (0, +00), in the perturbative RG analysis. ■^^•^^ However, 
with this method where the cosine and sine terms are treated 
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expansion of the critical values for small A/t: 







U/t 



FIG. 1 1 : Phase diagram of the half-filled extended Hubbard model at 
A/t = 0.1. The double line denotes the first-order transition, while 
the single lines denote the second-order transitions. 



perturb atively, we cannot see the correct picture of the DSG 
theory with the double-well potential structure which leads 
to the Ising transition. Instead, here we estimate the criti- 
cal value for the Ising transition from the semiclassical ar- 
guments: The critical value is determined from the condition 
GIIG\ = 1/4. 

We have used the above scheme to obtain the phase diagram 
shown in Fig.^] for which A/i = 0.1. The phase diagram 
at large U and V is similar to Fig. whereas a qualitative 
charge in the phase diagram is found in the region V^V "^t- 
In agreement with Fabrizio, Gogolin, and Nersesyan??^ we ob- 
tain two critical points ([/ci < C/c2) separating three phases on 
the V axis: the BI state, the BCDW state (= the SDI state^^), 
and the SDW state. From comparison of Figs. and ^2 we 
see that the BCDW state in Fig. has evolved continuously 
into the BCDW state when the alternating site potential A is 
switched on. The phase diagram in the IS.-V plane is shown in 
Fig-El where U jt = 1. Both A and V promote the BI state, 
while the SDW ground state is obtained for small A(<C C/) 
and 1^(< U). We find that the region of the BCDW state 
obtained in the EHM at A = is connected to the region of 
the BCDW state in the Hubbard model with alternating site 
potential at F = 0. 



Let us discuss in more detail the critical regime in the limit 
of small [/, V , and A. In this region we can safely neglect 
the irrelevant terms and set Gcs{l) = Gps{l) = in the RG 
equations (|^^}-(|^^}. First we consider the case V = 0. In- 
tegrating out the RG equations (|^^}-(|^^} analytically and 
following the criterion discussed above, we obtain asymptotic 



r7° 



\n{t/A) 

2TTt 

ln(i/A) 



C 



1 + C 



In(VA) 
,lnln(t/A) 



(4.21) 



ln(i/A) 



O 



ln(i/A) 



(4.22) 



where C and C are positive constants of order unity. 
The A dependence of U^-^ is different from the result in 
Refs. '32' since the lowest correction to 27ri/ln(i/A) is not 
0(ln[ln(VA)]/ln(</A)), but 0(l/ln(t/A)). Our results 
suggest that the ratio of U^2 C^ci becomes U^2/^ci = 
1 + C"ln[ln(t/A)]/ln(i/A). At present we do not know 
where this difference comes from. We extend this analysis 
to the case with finite U) and examine the V depen- 

dence of Uci and Uc2- We note that Gp{l) 7^ Gc{l) in this 
case since the SU(2) symmetry of the charge sector is broken. 
We integrate the RG equations analytically for small V ^ 
and obtain the corrections to order V, 



Uc2 = U^2 - V 



2 / 1 

3 ^ yn(VA) 

2 / lnlnft/A) \ 

3 ^ V ln(t/A) ; 



(4.23a) 
, (4.23b) 



implying that the BCDW state survives upon inclusion of the 
F(<C U) term. We note that Ud and Uc2 have a similar linear 
dependence on V. From Eqs. ( I4.23> and Figs.llllandll2l we 
conclude that the phase diagram exhibits reentrant behavior 
as V increases from zero with A and U being fixed at values 
near a quantum critical point. 

Since the Hamiltonian H' has three free parameters (U/t, 
V/t, and A/t) at half filling, the ground-state phase diagram 
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FIG. 12: Phase diagram of the half-filled extended ionic Hubbard 
model on the plane of A/t and V/t, where U/t = 1. 
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(a) 



(b) 





FIG. 13: Schematic phase diagram of the half-filled extended Hub- 
bard model at (a) 1/ = 0, (b) 1/ < t, and (c) 1/ > The single 
lines represent second-order transitions, and the double line in (c) 
represents a first-order transition. 



(a) 



BI 



BCDW 





FIG. 14: Schematic phase diagram of the half-filled extended Hub- 
bard model at (a) U <ti t and (b) U t. The single lines represent 
second-order transitions, and the double line represents a first-order 
transition. 



becomes a three-dimensional (3D) diagram. Instead of draw- 
ing such a 3D plot, here we show two-dimensional tomo- 
graphic phase diagrams. Figure ^] shows schematic phase 
diagrams in the A-C/ plane for three typical cases V/t = 0, 
V/t <^ 1, and V/t':^ 1. We see that the nearest-neighbor re- 
pulsion enhances the BI phase and destroys the BCDW phase 
at large V, where the direct transition between the BI and 
SDW phases is first order The recent numerical study of the 
ionic Hubbard model'"' reports a similar phase diagram as Fig. 
I13r a). The first-order transition line in Fig. lOf cl asymptoti- 
cally approaches the line U = 2A + 2V . 

Figure shows schematic phase diagrams in the IS.-V 
plane for U/t < 1 and U/t > 1. At large U and V there ap- 
pears a direct first-order transition between the BI and SDW 
phases in Fig. ll4r b). This first-order transition is in agreement 
with the results obtained from the strong-coupling analysis'*^ 
and numerical calcuIations42i^ 



C. Discussions on previous numerical results 

As mentioned in Introduction, many groups have already 
reported on numerical studies of the ground-state phase dia- 
gram of the ionic Hubbard model. Various numerical tech- 
niques were used in these studies, including the density- 
matrix renormalization-group (DMRG) methodpi^^i^^iSii^ 
the quantum Monte Carlo method r^''*' a finite-size cluster 
method,^** and a level crossing analysisiiS The main issue here 
is whether or not the SDI phase (BCDW phase) exists, and 
so far these numerical studies do not seem to have reached 



complete agreement yet. Although most of recent studies re- 
port that the SDI phase appears near the boundary between 
the SDW phase and the BI phasej?^'?^ii^A?'^i49i4i,4? there are still 
some conflicting claims in the literature. A less controversial 
issue^-' is the determination of the second critical value Uc2 at 
which a spin gap closes and which can be estimated by com- 
puting the spin gap directly^'*-^^ or by examining the BCDW 
order parameter. The determination of the critical point 
Uci and the critical behaviors around it are more controversial 
issues. One way to estimate the critical value Uci is to use 
the complex parameter introduced by Resta and Sorrela.Sft Its 
diverging behavior at [/ = Uci indeed allows one to deter- 
mine the critical point. Another way to determine the crit- 
ical point is to find a gap closing point in excitation spectra. 
Since the charge sector is responsible for the quantum phase 
transition ?aU — Uci, one might try to look at a charge gap 
directly. However, numerical studies have found that a naive 
charge gap does not vanish at the critical point and is always 
finite. Recent studies have shown^^'^^ "*^ that the excitation gap 
that vanishes at ?7 = Ud is the gap to the first excited state 
that has the same charge and spin quantum numbers as the 
ground state. Let us discuss this point in more detail below. 

In numerical studiespi^Sf^ the "charge gap" Ac was de- 
fined as Ac = Efi{L/2 + l,L/2) + £^0(^2 - 1,^/2) - 
2Eq{L/2,L/2), where Efi{N'^,Ni) is the lowest energy of 
a finite-size system with an even number of sites L that has 
up-spin and Ni down-spin electrons. This quantity Ac 
measures the energy of the excitation with the topological 
charge Q = ±1 and Sz — ±1/2 [Eq. ( 14. 7H . and is rather 
a single-electron excitation gap. According to the bosoniza- 
tion theory (Sec. IV A), the charge transition at ?7 = Ud is 
described by the "(^^" theory and is in the Ising universal- 
ity class. The transition occurs when two degenerate local 
minima of the effective potential for the charge fields merge 
into a single local minimum. As one approaches the transition 
point from the Ising ordered phase (that is, the SDI phase), 
the topological charge Q — ±2ao/n of a lowest-energy exci- 
tation is decreasing to zero, while excitations with Q = \ re- 
main massive. Therefore the charge gap Ac does not vanish at 
this Ising critical point, and this quantum phase transition can- 
not be detected with Ac. Qin et al. and Manmana et al. also 
used Ac = El (L/2, L/2) - Eq (L/2, L/2) in their numeri- 
cal analysis, where Ei (7V| ,Ni) is the energy of the first ex- 
cited state. •'^•'^^ The quantity Ae measures excited states with 
the same number of electrons, whose total topological charge 
Q = in the sine-Gordon scheme. In the Ising ordered phase, 
the first excited state with the topological charge Q = would 
be a bound state (or breather) of a soliton with the topological 
charge +2ag/TT and an antisoliton with the charge — 2ae/7r, 
whose energy vanishes at the critical point. On the other hand, 
in the Ising disordered phase near the critical point, the po- 
tential is almost flat and has very small curvature. The low- 
energy excitations would then be small oscillations around po- 
tential minima (rather than soliton/antisoliton) whose energy 
approaches zero as U ^ Ud — 0. Thus the exciton gap Ac 
is a right measure to detect the quantum phase transition at 
U = Ud- 



16 



TABLE III: Possible ground states and the position of locked phase 
fields, determined from Eq. <5.4> . 



Phase 



SDW 
CDW 

PI (for gs > 0) 
PI (for gs < 0) 
BSDW 



(O,±70),(7r,±(7r-7^)) 
(±7,,0),(±(7r-79),7r) 

(0,0),(7r,7r) 

(0,7r),(7r,0) 
(±V2,±^/2) 



V. EFFECT OF BOND DIMERIZATION 

In this section, we consider the ID EHM with staggered 
bond dimerizationi^i^ i.e., the Peierls modulation of the hop- 
ping matrix element. The total Hamiltonian H" is given by 
H" = H + Hs, where H is defined in Eq. (ED and 
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FIG. 15: Positions of locked phase fields i 
states for qs > 0. 



and (j) in the respective 



H.c. 



(5.1) 



Without loss of generality we can assume S > 0. When V = 
0, the model is called "Peierls-Hubbard model." The one- 
dimensional Mott insulator, realized when U > and = 0, 
is known to be unstable against the Peierls distortion,'-^ and 
as a result the ground state changes from the SDW state into 
the BCDW state regardless of the magnitude of the Hubbard 
interaction U. Such an instability comes from the fact that 
the bond dimerization tends to concentrate the electron den- 
sity onto bonds, without any conflict with the Hubbard, U, 
repulsion.'"' However, the nearest-neighbor Coulomb repul- 
sion V competes with this S term, since the V interaction 
likes to localize two electrons on a single site and promotes the 
CDW state. Here we investigate the instability of the BCDW 
state against the intersite Coulomb repulsion V, and clarify 
the critical behavior near the transition between the BCDW 
state and the CDW state. 

The bond dimerization Hs is bosonized as Hs — J dxHs, 
where 



A. Semiclassical analysis 

We begin with semiclassical analysis of the model with the 
gs term. We neglect spatial variations of the phase fields in 
H + Hs and consider the potential 

- gs cos 2(j) — gcs cos 29 cos 20 
•r.«^ (5.3) 



Vs{eA) = 



-(?c cos 29 - 
" gs cos 9 cos < 



where = 53|| < 0. 

First, we consider the simpler case where gcs — 0, which 
corresponds to the situation where gcs is irrelevant in the RG 
sense. The potential in this case is 



-gc cos 29 + gs cos 2(f) — gs cos 9 cos 0. (5.4) 



The positions of the potential minima are determined 
by the saddle-point equations, dVs {9 , (j)) / 09 = and 
dVs{9, = 0. We find that the potential has the double- 

well structure for the 6 (0) phase field when gc < —\gs\/4: 
(gs > Iga 1/4). Here we introduce 7^ and 7^ (0 < 7^, 7^ < tt) 
defined by 



Ti-s 



gs 



2(7ra)- 



■ cos f cos ( 



(5.2) 



and gs — SirSa. One finds that the EHM with the bond dimer- 
ization also has a two-component DSG structure. Here the 
charge phase field 9 is subjected to the potential cos 9 instead 
of sin 9 of the term [Eq. ( 14. 2H . while the locking potential 
for the spin phase field has the same structure as that of the 
gA term. 

It is important to note that the BCDW order parameter 
C'bcdw takes a nonvanishing expectation value for any U and 
y if (5 ^ 0, as Hs oc Obcdw- In this section we will not use 
the term BCDW to characterize phases, and, in particular, the 
phase containing the trivial Peierls insulator (U — V = and 
(5 ^ 0) is called the Peierls insulating (PI) phase. 



7.° 



9s 



4 



gs 
'^gs 



(5.5) 



for \g5/gc\ < 4 and \gs/gc\ < 4, respectively. The solutions 
to the saddle-point equations can be classified into the follow- 
ing four classes: (i) the PI state, [9, (j)) = (0, 0), (0, tt), (tt, 0), 
or (tt, tt) [for gs > 0, the phase fields are locked at {9, (f) = 
(0, 0) or (tt, tt), while for 55 < the phase fields are locked as 
(6*, (j)) = (0, tt) or (tt, 0)]; (ii) the pure BSDW state, (6*, 0) = 
(7r/2,±7r/2) or (-7r/2, ±7r/2); (iii) the "SDW" state with 
both the SDW order and the BCDW order, (6*, 0) = (0, ±7^ 
or (tt, ±(7r - 7^)); and (iv) finally, the "CDW" state with 
both the CDW order and the BCDW order, {9, <j)) = (±7^, 0) 
or (±(7r — 7g), tt). The possible ground states and positions 
of locked phase fields are summarized in Tablellllland Fig. 1151 
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FIG. 16: Phase diagram obtained by minimizing the potential energy 
Vs{9, 0) [Eq. OZJ]- The phase boundary of the BSDW state is given 
by the curve QcQs = — (/|/16 with Qc < 0. The phase boundary 
between the PI state and the SDW state and that between the PI state 
and the CDW state are given by the lines Qs = |5i|/4 with > 
|(7a|/4 and gc — ~\gs\/4: with gs < lff«|/4, respectively. All the 
phase transitions in this figure are continuous. A multicritical point 
is at (gc,gs) = {-\gs\/4.,\gs\/4). 



In these states the potential energy reads 

V° - 



BSDW 



V" 

''CDW 



9c + gs- \gs\, 
9c - 9s, 



~9c - 9s 



9l 



9} 

'9c + 9s + -§^- 



(5.6a) 
(5.6b) 

(5.6c) 
(5.6d) 



In deriving Eqs. (I5.6c> and (I5.6d> . we have assumed \gs/gs\ < 
4 and \gs/gc\ < 4, respectively. The PI state is stabilized 
by the first-order contribution of the gs term. Furthermore, 
if 5s > (gc < 0), the SDW state (the CDW state) is also 
stabilized due to second-order contribution of gs. The phase 
diagram obtained by comparing these energies is shown in 
Fig. EH 

From the above semiclassical analysis one might conclude 
that the topological charge Sz [Eq. ( 14. 7H becomes fractional 
in the SDW phase and that the Ising-type phase transition in 
the spin sector takes place on the boundary between the PI 
state and the SDW state. However, as discussed in Sec. IIVI 
the global SU(2) symmetry prohibits the Ising criticality in 
the spin sector and changes the SDW phase in Fig.[^into the 
PI phase. 

Next we include the gcs term. Table |III] still stands if we 
replace gc and g^ with gc - \gcs\ and gs + \gcs \ in 7g and 7^, 
respectively. The phase diagram obtained by minimizing the 
potential energy Vs{d, </>) is shown in Fig.^] New features 
compared with Fig.^]are the appearance of a first-order tran- 
sition line and of the new phase in which the ground state has 
the coexisting order of the SDW, CDW, BCDW, and BSDW. 



BSDW 



SDW 




FIG. 17: Phase diagram obtained by minimizing the potential en- 
ergy Vsie,4>) [Eq. O)] drawn for \gs\/8 < \gcs\ < |3i|/4. 
Multicritical points are located at (ffcffs) = {~\gcs\,\gcs\) and 
(-ils^i + Iffcsl, j\gs\ - \9cs\)- The boundary of the BSDW phase 
is igc + \gca\){gs — \gcs\) < —gl/lQ. The edges of the PI phase are 
defined by the lines = -ils^l + lffcsl andg^ = ilStfl-lffcsj. The 
double line denotes a first-order transition, and the single lines denote 
continuous transitions. Within the semiclassical analysis the ground 
state in the shaded region has the coexisting order of the SDW, CDW, 
BCDW, and BSDW. 



The new phase is shown as the shaded region in Fig.[n] which 
is surrounded by the three curves defined by 



^2 

{9c + \9cs\){gs - \gcs\) = 

16 

5| 



{9c + \9cs\){gs + \gcs\f 



16 

9l 



(5.7a) 



(5. - \9cs\), (5.7b) 



{9s - \9cs\){gc - \gcs\Y = -^{9c + \9cs\). (5.7c) 

Let us focus on the phases which can be realized when gs — 
gc, in view of the fact that in the extended Hubbard model 
both gs{^ 9i±) and gc{= g3±) are given by {U - 2V) in 
the lowest order Along the line gs ~ gc in Figs. 1161 andll7l 
there are three possible phases: the SDW state, the PI state, 
and the CDW state. Since the SDW state is prohibited by the 
SU(2) symmetry and becomes the PI state, we expect to have 
only two phases, the PI state and the CDW state, and a single 
phase transition between them. The transition is continuous 
at \gcs/9s\ ^ 1 and changes into a discontinuous transition 
when gcs exceeds |55|/4. 



B. Renormalization-group analysis 

Next we perform perturbative RG analysis to take into ac- 
count quantum fluctuations. The one-loop RG equations for 
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coupling constants in 7i + Tis are given by 

— = + Gs + -Gs Gp + Gs Gc 

3 3 3 

— - Gs Gs + - Gs Gcs ~ g ^p*' (^-^^ 

— Gp = + — G^ + 2 G^ + G^g + Gs Gps, (5.9) 

^Gc = + — G^ + 2 Gp Gc — Gs Gcs ~ Gcs Gps, 

(5.10) 

— Gs = — — G^ — 2 G^ — Gc Gcs ~ G^j,, (5.11) 

^Gcs = +\Gl-2Gcs + 2GpGcs~^GsGcs 

2 Gc Gs 2 Gc Gps 4 Gcs Gps, (5.12) 

^Gps = — — G^ — 2 Gps + 2 Gp Gs 

-4GcGcs-4GL-4GsGps. (5.13) 

The initial value of Gs{l) is given by Gs{0) — 26/t and 
those of the other coupling constants are Gi,{0) — g^/ (47rt). 
We note that these RG equations are invariant under the sign 
change of Gs{l)- We can thus assume Gs{0) > without 
losing generality. 

To find the ground-state phase diagram of the system, we 
solve the scaling equations ( l5.8> -( l5T3t numerically, as in the 
preceding sections. We determine to which phase the ground 
state belongs by looking at which one of the couplings Gs{l), 
Gc{l), Gs(/), and Gcs{l) becomes most relevant. For repul- 
sive U and V there are four possibilities as listed below. 

(i) If Gc is most relevant and Gc{l) 1 at Z = Zp+, then 
the phase field is locked at 9 ~ or vr, and the effective 
Hamiltonian for the spin sector at I > lp+ becomes 

= '^Y.^dAvf -^-^G:{dA+){dA-) 
Zn ^ — ' TT 
P=± 

T^G| cos0+-^G: COS20, (5.14) 

where G* = Gs{lp+) - Gcs{lp+) and G| = Gsilp+), and 
the sign — /+ of the G| term corresponds to the location of 
the phase locking 9 = O/tt. This effective theory is the same 
as Eq. (I4.19> . As seen before, regardless of the sign of G*, 
the phase field (p is locked at = or vr depending on the 
position of the charge phase locking = or vr. Thus we 
have the phase locking {6,(f)) — (0,0) or (7r,7r), i.e., the PI 
state as the ground state. We note that due to the SU(2) spin 
rotation symmetry the SDW state cannot be realized even if 

g; > 0. 

(ii) If Gc is most relevant and Gc{l) — 1 at Z = Zp_, 
then the phase field 6 is locked at 6 = ±tt/2. The effective 
Hamiltonian for the spin part is 

^f- = ^Y.'^d,'^p^'-^G:{d,cj^+){d,q^^) 

ZTT ^ — ' TT 

p 

+ — ^G!cos20, (5.15) 



where G* = Gs{lp-) + Gcs{lp-)- We have verified numer- 
ically that G* always becomes negative in this case. The 
G* terms are then marginally relevant [Gl{l) — > —oo un- 
der scaling]. The phase fields are then locked at {9,(p) = 
(±7r/2, 0), (±7r/2, tt), which corresponds to the CDW phase 
with 76/ tt/2 (i.e., gs 0, see Table HlH . Since Hs oc 
Obcdw, the order parameter of the BCDW should have a 
nonvanishing expectation value. We thus conclude that the 
ground state is in the CDW phase. 

(iii) If either Gs or Gcs is most relevant, both charge and 
spin fluctuations are suppressed. In this case the semiclassical 
treatment is justified, and we can determine to which phase 
the ground state belongs by substituting Gc and Gs to gc and 
gs in Fig. [HI 

(iv) If Gs is most relevant and Gs{l) — 1 at / — l^, 
the spin fluctuations are suppressed and the phase field cj) is 
locked at ^ or tt below this energy scale. The effective 
Hamiltonian at I > Ic is given by 

ZTT ^ — ■ TT 

p=± 

T ^ G*s cos9- ^G*c cos29, (5.16) 
TTa^ TTa^ 

where G; = Gp{Q - Gps{la), G* = Gc{la) + Gcs{L), 
and G| = Gs{la)- The sign — /+ of the GJ term corresponds 
to the phase locking (p = O/vr. Both of the nonlinear terms 
cos 9 and cos 29 are relevant perturbations. If G* < 0, these 
two terms compete with each other, and this DSG model ex- 
hibits the Ising criticality. The ground state is either in the PI 
phase or in the CDW phase, and there is an Ising-type quan- 
tum phase transition between the two phases. Here we esti- 
mate the Ising critical point from the semiclassical analysis. 
That is, the critical value is determined from the condition 
G*/G*s = -1/4 (see Fig.O. If G* > 0, these two terms 
do not compete and thus the phase locking is 6* ~ {tt) for 
(j> = (tt), where the ground state is the PI state. 

The resultant phase diagram in the U-V plane is shown in 
Fig.^] In the weak-coupling region, the transition from the 
PI state to the CDW state is characterized by the appearance 
of the double-well structure of the effective potential to the 
9 field, and thus the phase transition in Fig. ^] belongs to 
the Ising universality class. As we increase U and V, there 
appears a tricritical point at {Uc, Vc) ~ (4.9t, 2.3t), where the 
phase transition changes from second order to first order 

Figure[l9]shows schematic phase diagrams in the S-U plane 
for V and V/t ^ 1. When 6 — 0,we obtain three phases 
(the CDW, BCDW, and SDW phases) for < i (a) and two 
phases (the CDW and SDW phases) for y > i (b), as we 
discussed in Sec.|ni](see Fig.0i. Upon turning on S, the SDW 
ground state changes into the PI state, where the transition 
is described by the Gaussian theory. On the other hand, the 
BCDW state changes into the PI state without accompanying 
any singularity: This change is merely lifting of the doubly 
degenerate BCDW ground states. 

Figure|20|shows schematic phase diagrams in the S-V plane 
for [7 = 0, U/t < 1, and U/t » 1. At [/ = we have 
a single critical value Vc which has the S dependence given 
by Vc oc l/ln(t/(5) for small 6. As U and V increase, the 
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(a) 



(b) 



(c) 



FIG. 18: Phase diagram of the half-filled extended Hubbard model 
with ijt — 0.1. The second-order transition line (single line) turns 
into the first-order transition line (double line) at the tricritical point 
{Uc,Vc) ^ (4.9t,2.3t). 



phase boundary approaches the U — 2V line. The asymp- 
totic form of Vc for U,V ^ 6 and S t is given by 
V; = i[/ + C" U{5/tf^l''\ where C" is a numerical con- 
stant of the order of unity (see also Fig.FTSt. 



VI. CONCLUSIONS 

In this paper we have studied the ground-state phase dia- 
gram of the one-dimensional extended Hubbard model with 
on-site and nearest-neighbor repulsion U and V . By includ- 
ing higher-order corrections to coupling constants in the g- 
ology, we have given a plausible theoretical argument within 
the RG approach for the mechanism of the appearance of the 
BCDW phase at C/ sa 2V in the weak-coupling limit. Our 
two-step RG approach, however, is not complete in that there 
remains a weak cutoff dependence in the phase boundaries. 
This, albeit minor, defect should be resolved with use of a 



(a) 



~SDW 



PI 





FIG. 19: Schematic phase diagram of the half-filled extended Hub- 
bard model at (si)V <^ t and (b) V ^ t. The single lines represent 
second-order transitions, and the double line represents a first-order 
transition. 



'-•BCDW 
-SDW 



FIG. 20: Schematic phase diagram of the half-filled extended Hub- 
bard model at (a) [/ = 0, (b) [/ < t, and (c) U '> t. The single 
lines represent second-order transitions, and the double line repre- 
sents a first-order transition. 



more sophisticated systematic RG procedure. Away from 
the weak-coupling limit the umklapp scattering between the 
parallel-spin electrons 53|| tends to destabilize the BCDW 
state and eventually gives rise to a bicritical point where the 
two continuous-transition lines merge into the SDW-CDW 
first-order transition line (Fig. 0. We should note, however, 
that there still remains a difficult question as to whether our 
phase diagram is qualitatively correct near the multicritical 
point (which we call bicritical). One could imagine, for exam- 
ple, a possibility that a continuous phase transition between 
the BCDW state and the CDW state becomes first order be- 
fore reaching the multicritical point, due to higher-order ef- 
fects that are ignored in our analysis. If the coiTect topology 
of the phase diagram is indeed the same as ours (Fig.Q^ then 
the critical properties of the multicritical point remain to be 
understood. We hope that these issues will be resolved by fu- 
ture studies. 

We have also examined effects of additional staggered site 
potential and bond dimerization in the extended Hubbard 
model. In the presence of the staggered site potential, we have 
found that the BCDW state is smoothly connected to the SDI 
phase which is obtained for V — Qhy Fabrizio et al.^^ In this 
BCDW phase the BCDW order coexists with the CDW order, 
and the quantum phase transition between the BI phase (or the 
CDW phase) and the BCDW phase belongs to the Ising uni- 
versality class (c = i CFT). For finite V the BCDW phase 
is also destabilized by the 93 1| term, and the direct first-order 
quantum phase transition between the SDW state (= Mott in- 
sulating state) and the BI state takes place (Fig. II 1> . In the 
presence of the staggered bond dimerization the SDW phase 
becomes unstable and the ground state at V" = turns out to 
be the Peierls insulating state. For V ^ the phase diagram 
consists of two phases, the PI state and the CDW state, which 
are separated by a phase transition line of the Ising criticality 

(Fig.d}. 
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APPENDIX A: BOSONIZATION 

In this section, we derive the phase Hamiltonian of the ID 
extended Hubbard model by using the Abelian bosonization 
method.''^ We include not only the marginal terms but the 
leading irrelevant terms which play a crucial role in the first- 
order SDW-CDW transition at strong coupling. 

The Lagrangian for the free massless boson theory in a two- 
dimensional Euclidean space is given by 



Lb 



= i dx 



V 



(Al) 



where 6* is a bosonic field, r is the imaginary time, and v is 
velocity. The variable canonically conjugate to 9 is given by 



n 



dL 

'd9 



2'KV 



dr9, 



(A2) 



where 9 = d9/dt and t is the real time (t = it). As usual this 
system is quantized by imposing the commutation relation at 
equal times: [9{x),'n{x')] — iS{x—x'). Thus the Hamiltonian 
for the free boson theory is given by Ho — i J dxIldT9 + Lg, 



I.e., 



— I dx 

An 



(A3) 



Introducing two copies of this theory with fields 9 and (j) 
and velocity v = vp, we arrive at Hq [Eq. ( I2.11H . where the 
fields 9 and represent the "charge" and "spin" degrees of 
freedom. The chiral bosonic fields 9± (a;, t) and 0± [x, r) are 
introduced in Eqs. (|^^) and ( I2.13> , respectively, where the 
right-moving (left-moving) fields are functions of T — i{x/vF) 
[r + i{x/vF)]-^^ The phase field 9 (0) and its dual phase field 
9 {(p) are written in terms of the chiral fields as 



6l_, 9 



(A4) 
(A5) 



They satisfy the following commutation relations: 

[9{x),9{x')] = [(j){x),4>{x')] = -i2TTe{-x + x'), (A6) 

where Q{x) is the Heaviside step function. 

The electron field operators 1/;^ ^ (x) are given in Eq. j2.17> 
in terms of a new set of chiral bosonic fields ipp,a- introduced 
in Eq. ( I2.15> . In this bosonization scheme ip+.a- and V'-.o- an- 
ticommute, and we only need to introduce the Klein factor 
Ka to ensure the anticommutation relation between fields with 
different spins; cf. the so-called constructive bosonization 
method.''^ From Eqs. ( I2.16t and ( I2.17> the electron-density 
operator becomes 



Pp, 



,{x)= : ^pl^^ipp^„: 



1 d 
27r dx 



(A7) 



As is well known, the Hamiltonian density of free bosons 



VF 

An 



dx 



p,CT 



is equivalent to the Hamiltonian density of free fermions with 
linear energy dispersion, Eq. ( 12. 9> . This can be shown, for 
example, by using the OPE method. 

Next we bosonize the interaction term Hint- Without the 
nearest-neighbor repulsion V, this can be easily done asJ^j^ 



V=0 _ 94\\ + 94± 



47r2 
94\\ - g4_L 
47r2 

02 II + 92± 



{d^9+f + {d,9^f 



.9111 



27r2 

5211 - 52± - 511 



27r2 



53-L „^ , .91 

COS 29 + 



2{-Ka) 



{dj+){dj^) 

{dx(j)+){dx(p-) 

cos 2(f>, (A9) 



2(7ra)2 



where g's are given in and below Eq. ( 12. St . In the presence 
of V, the matrix element of the umklapp process with parallel 
spins iJg^i^ [the g^ process in Eq. ( 12. 4H has a finite amplitude 
at lowest order in g-ology. This term can be bosonized as 



^3311 — 



53|t 



2(7ra)- 



■ cos 29 cos 20, 



(AlO) 



where g^ = —2Va in the lowest order in V. This term, which 
couples the charge and spin degrees of freedom, is often ne- 
glected since it is an irrelevant perturbation with scaling di- 
mension 4, consisting of dim[cos 29] = 2 plus dim[cos 20] — 
2. Cannon and Fradkin were the first to suggest that this term 
should play an important role in the first-order SDW-CDW 
transition in the half-filled EHM.^^ Voit then derived RG equa- 
tions including this term. However he did not include all the 
operators with scaling dimension 4 and failed to keep the spin- 
rotational SU(2) symmetry.^'* We have to be careful in deal- 
ing with the V interaction to include the important terms with 
scaling dimension up to 4. To this end, we focus on the V 
interaction and bosonize each scattering process separately. 

First, the term^ representing the backward scattering 
with parallel spins is bosonized by using Eq. il.lli as 

Va ^ 0t ,^ (x) ip-p^a (x) i^lp^^ {x + a)il;p^a{x + a) 



Va 



(27ra)2 



^ip[9(x^a)—0{x)]-^ips[(l)(x-\-a)—<p{x)] 



Va 
2^ 



(d.Opf + 2 [3^9+) (9, 



Va 
2^ 



^a^ (9,0)2(9^0)2 + 



(All) 
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where we have expanded the exponent in the second line 
up to the order a* for the 9 sector and the (f> sector, sep- 
arately. Since we are interested in operators that couple 9 
and (j> as in Eq. ( lAlOt . we have discarded dimension-4 terms 
such as a*{dxd)'^ and a*{dx4')* that involve only one sec- 
tor. Such terms as {dxd+){dx9-) and {dx<j>+){dx4i-) are al- 
ready retained in Eq. ( IA9l l. while the last term proportional 
to {OxfiY {dxi'Y is ^ new term with scaling dimension 2 + 2, 
which was missed in Ref. |3 We note that the Fermi ve- 
locity is renormalized by the gi\\ term due to the presence of 

Tlipi'^x^ipY ^iid Tlipi'^x't'pf' ■ This is in contrast with the con- 
ventional treatment where the velocity renormalization comes 
only from the forward scattering term g^w^ 



In a similar way, the interaction terms of backward and 
umklapp scattering with opposite spins (so-called gi±_ and 
termspi respectively) are bosonized as 



^ip[O{x+a)~0{x)]-ips[4>{x+a)+^{x)] 



Va 
" (27ra)2 

.J}^cos2(f, + ?^{dx9+){dx9-) cos 20 



D,S — ± 



2Va 



COS 20 + ■ 



(A12) 



'^a^il)l „{x) il^-p.aix) i/jI^x + a) llj-p,a{x + a) 



E 



'iple{x+a)+e{x)]+ipsl<t)(x+a}~4>{x)] 



Va 

(2^ , 

^ 00329 — {dx(j)+) {dx(j)-) cos29 



2(7ra)2 

2Va 
~ 4^ 



27r2 



E(^^ 



xHJpj 



COS 20 H , 



(A13) 



Hence the total Hamiltonian is given by 



H = ^ E [^p {dx9pf + v„ {dx<f)pf 



p 



cos 20 + „ „ ^ COS 20 



27r2a2 
Fa 



27r2a2 
cos 20 cos 20 



+ 



Va 
2^ 

Va 
2^ 



Va^ 
4^ 



Y,{dx9pf + 2{dx9+)(dx9_) 
p 

Y^idxCl^pf +2{dx<P+){dx<j)-) 

p 

Y^idxOpf + 2{dx9+)idx9-) 
- p 

E(5.0p)' + 2(a,0+)(a,0_) 



cos 20 



cos 20 



(A14) 



The renormalized velocities are given by Vp = 2ta + {U + 
6V)a/{2TT) and v^r = 2ta - (U - 2V)a/{2TT). The coupling 
constants gi± and g3± are defined in Eq. (12. 8> . and .gp(= g2\\ + 
92± - 9i\\) andg„(= -5211 + 92± + 9i\\) are given by 



gn = (U + 6V)a+ —(U -2Vfa+—V'^a, (A15a) 
47ri 7rt 

g^ = (C/-2y)a- -^([/-2y)2a- — (A15b) 
47rt Tri 



For the discussion of the SDW-CDW transition in the ID 
EHM, it is sufficient to have the coupling constants of di- 
mension 4 in lowest order in V. We note that due to the 
SU(2) spin-rotation symmetry of the theory, the coupling con- 
stants for spin degrees of freedom must satisfy g^r — gi±, 
in any order of U and V. To proceed further, we neglect 
the te rms that involve V Y.pi^xOp)'^ or V J2pi9x4'p)^ in Eq. 
JA14I I. These terms can lead to renormalization of the velocity 
through the RG transformation (see Appendix|nji. This effect 
can be ignored if we are interested in qualitative feature of the 
ground-state phase diagram of the model. The final form of 
the bosonized Hamiltonian is thus given by Eq. j2.19> . 



where a =| (t) for a =t (i). The potential cos 20 in Eq. 
( IA12> and the potential cos 20 in Eq. ( IA13> are akeady re- 
tained in Eq. (IA9> . while the other terms are new and have the 
scaling dimension 2 + 2. 



The forward-scattering terms (g2\\, .g2±> g4\\, and 54^) do 
not generate operators of dimension 2 + 2. 



APPENDIX B: DERIVATION OF 
RENORMALIZATION-GROUP EQUATIONS 

In this section, we derive one-loop RG equations for the 
coupling constants including those operators with higher scal- 
ing dimension. Our derivation is based on the operator prod- 
uct expansion (OPE) method. The interaction part of the ac- 
tion 5*1 in the presence of the staggered site potential Ga is 
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given by 



Si 



^ / d\ {d,B) {d,e) I d\ (9,0) 



Ga [ d?r 



: sm y : : cos o : 



: cos 26 : 



: cos 2<p : 



TT a 



: cos 20: :cos2(/): 



J d^r {dj) {d,e) ■.cos2(j): 
+ ^ /"d^r :cos26': 



+ ^ / dVa^ (9,0) (9,0) (9,</.) (9,(/.) , (Bl) 



where z = vpT + ix,'z = vpT — ix, d'^r = vpdxdr, and 
Gi = gi/2'KVp. In this section, the operators are exphcitly 
normal ordered. 

In order to derive the RG equations, we use the following 
OPE's: 



Jp{z) Jp{w) = 
Jp{z) Jp{w) 



1 



{z — w)-^ 

1 



{z — wy 



Jp{z) :e 
Jp{z) :e' 



2{z - w) 

—a 
2{z-w) 



(B2a) 
(B2b) 
, (B2c) 
, (B2d) 



2\z\ 
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where we have introduced U(l) currents: Jp{z) = idzO(z, z), 
Jp{z) = -idzO{z,z), J<j(z) = idz(t){z,z), and Ja{z) = 
—idz(j){z, z). The parameters a and /3 (a + /3 ^ 0) in the 
vertex operator are numerical constants which determine the 
scaling dimension. In deriving the above OPE's, we have 
used the Wick theorem and the correlators: {0+{z) 0+(w)) — 
-iln(z - w), (0_(z)0_(w)) = -iln(z - w), and 
{d{z,z) 9{u!,0j)) — — ln|z — From Eg. (IB2> . one finds 
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:cos[(a + /3)0]: 
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Exchanging 9 ^ (p and p ^ a yields the OPE's for spin phase 
fields. 



Expanding the action in powers of coupling constants and 
integrating out short-distance parts, we obtain the scaling 
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equations, 
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Here we note that the number of the RG equations can be 
reduced due to the spin-rotational SU(2) symmetry. To show 
this point more transparently, we introduce X{1), Y{1), and 
Z{1) by X{1) = Gail) - Gs(0, Y{1) = Gcs(0 - Gc.(/), and 
Z{1) — Gps{l) — Gpa{l). Their RG equations are obtained 
from Eqs. (IE51- (IbT21 as 

= 2 Gs X + (Gc - Gcs) (B13a) 

dl 

j^Y = 2(-l + Gp + Gs + Gps)Y 

+ 2iGc-Gcs)iX - Z), (B13b) 

j^Z = -2(l-Gs)Z-2(Gp + Gps)X 

~'i{Gc~Gcs)Y. (B13c) 

One immediately finds that, if X(0) = Y{0) = Z{0) = 0, 
they vanish for all /, i.e., X{1) = Y{1) ^ Zil) = 0. 
This implies that G^(0 = Gs(0, Gcs(0 = Gca{l). and 
Gps (0 = Gpa (l), which are nothing but the constraints on the 
coupling constants due to the spin-rotational SU(2) symmetry. 
In this case, we can set Go- (0 — Gs{l), Gca (0 — Gcs (0 ' and 
Gpail) = Gps(0 in the RG equations ^4\-^T2\ . Then the 
RG equations are given by Eqs. ( I4.12> -( |4.17> . The RG equa- 
tions for the ID EHM without the staggered site potential are 
obtained by setting Ga(0 = 0, Eqs. ( I^0b -( ir24b . 

The RG equations can also be obtained in the presence of 
the bond dimerization in a similar way. 
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